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We present an analysis of exclusive charmless semileptonic B-meson decays based on 377 million 
BE pairs recorded with the BABAR detector at the T{AS) resonance. We select four event samples 
corresponding to the decay modes i?" — >■ tv' t^v, B"*" — ^ n'^l^u, B" — ^ p~£'^v, and B^ — >■ p'^l^u, and 
find the measured branching fractions to be consistent with isospin symmetry. Assuming isospin 
symmetry, we combine the two B — >■ ttEu samples, and similarly the two B — >■ plv samples, and 
measure the branching fractions BiB" -^ 7r"^+i/) = (1.41 ± 0.05 ± 0.07) x 10"'' and Z3(B° -s- 
p~ t^v) — (1.75±0.15±0.27) X 10"*, where the errors are statistical and systematic. We compare the 
measured distribution in q^ , the momentum transfer squared, with predictions for the form factors 
from QCD calculations and determine the CKM matrix element \Vub\- Based on the measured 
partial branching fraction for B — ^ ■niv in the range (f < 12 GeV'^ and the most recent LCSR 
calculations we obtain \Vub\ ~ (3.78 ± 0.131o'4q) x 10~^, where the errors refer to the experimental 
and theoretical uncertainties. From a simultaneous fit to the data over the full (f range and the 
FNAL/MILC lattice QCD results, we obtain \Vub\ = (2.95 ±0.31) x 10"^ from B -> ivlv, where the 



error is the combined experimental and theoretical uncertainty. 
PACS numbers: 13.20. He, 12.15.Hli, 12.38. Qk, 14.40.Nd 
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I. INTRODUCTION 

The elements of the Cabibbo-Kobayashi-Maskawa 
(CKM) quark-mixing matrix are fundamental parame- 
ters of the Standard Model (SM) of electroweak inter- 
actions. With the increasingly precise measurements of 
decay-time-dependent CP asymmetries in i?-meson de- 
cays, in particular sin(2/3) [l|,y], improved measurements 
of the magnitude of Vub and Vcb will allow for more strin- 
gent experimental tests of the SM mechanism for CP vio- 
lation [3] . This is best illustrated in terms of the unitarity 
triangle, the graphical representation of one of the unitar- 
ity conditions for the CKM matrix, for which the length 
of the side that is opposite to the angle /3 is proportional 
to the ratio |X^b|/|Kb|- The best method to determine 
\Vub\ and I Vcb I is to measure semileptonic decay rates for 
B — >■ Xcliy and B — )■ X^ii^ {Xc and X„ refer to hadronic 
states with or without charm) , which are proportional to 
|Vc6p and iKbP, respectively. 

There are two methods to extract these two CKM el- 
ements from B decays, one based on inclusive and the 
other on exclusive semileptonic decays. Exclusive de- 
cays offer better kinematic constraints and thus more 
effective background suppression than inclusive decays, 
but the lower branching fractions result in lower event 
yields. Since the experimental and theoretical techniques 
for these two approaches are different and largely inde- 
pendent, they can provide important cross checks of our 
understanding of the theory and the measurements. An 
overview of the determination of \Vub\ using both ap- 
proaches can be found in a recent review [J] . 

In this paper, we present a study of four exclusive 
charmless semileptonic decay modes, B^ — >• Tr~£~^i^, 
B+ -^ TT°i+v, B° -^ P'i^i^, and B+ -^ p°e+iy Q, and 
a determination of \Vub\- Here £ refers to a charged lep- 
ton, cither e"*" or fi'^ , and v refers to a neutrino, cither 
i^e or v^ . This analysis represents an update of an earlier 
measurement [y] that was based on a significantly smaller 
data set. For the current analysis, the signal yields and 
background suppression have been improved and the sys- 
tematic uncertainties have been reduced through the use 
of improved reconstruction and signal extraction meth- 
ods, combined with more detailed background studies. 
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The principal experimental challenge is the separation 
of the B — >■ Xuii^ from the dominant B — >■ Xcii^ decays, 
for which the inclusive branching fraction is a factor of 50 
larger. Furthermore, the isolation of individual exclusive 
charmless decays from all other B — ?► X„£i/ decays is diffi- 
cult, because the exclusive branching ratios are typically 
only 10% oiB{B -^ X^iv) = (2. 29 ±0.34) x lO'^ Q, the 
inclusive branching fraction for charmless semileptonic B 
decays. 



The reconstruction of signal decays in e+e^ — > 
T(4S') — ?► BB events requires the identification of three 
types of particles, the hadronic state X^ producing one or 
two charged and/or neutral final state pious, the charged 
lepton, and the neutrino. The presence of the neutrino is 
inferred from the missing momentum and energy in the 
whole event. 



The event yields for each of the four signal decay modes 
are extracted from a binned maximum-likelihood fit to 
the three-dimensional distributions of the variables toes, 
the energy-substituted i3-meson mass, AE, the difference 
between the reconstructed and the expected i3-meson en- 
ergy, and q^ , the momentum transfer squared from the 
B meson to the final-state hadron. The measured differ- 
ential decay rates in combination with recent form-factor 
calculations are used to determine |Vuh|. By measuring 
both B — >■ niv and B — )• piv decays simultaneously, we 
reduce the sensitivity to the cross feed between these two 
decay modes and some of the background contributions. 



The most promising decay mode for a precise determi- 
nation of \Vub\, both experimentally and theoretically, is 
the B — > ■kH.u decay for which a number of measurements 
exist. The first measurement of this type was performed 
by the CLEO Collaboration [3|. In addition to the ear- 
lier BaBAR measurement mentioned above [6| , there is a 
more recent BaBAR measurement [9[ in which somewhat 
looser criteria on the neutrino selection were applied, re- 
sulting in a larger signal sample but also substantially 
higher backgrounds. These analyses also rely on the 
measurement of the missing energy and momentum of 
the whole event to reconstruct the neutrino, without ex- 
plicitly reconstructing the second B-meson decay in the 
event, but are based on smaller data sets than the one 
presented here. Recently a number of measurements of 
both B -^ Txlv and B -^ piv decays were published, 
in which the BB events were tagged by a fully recon- 
structed hadronic or semileptonic decay of the second B 
meson in the event [10, llll • These analyses have led to a 
simpler and more precise reconstruction of the neutrino 
and very low backgrounds. However, this is achieved at 
the expense of much smaller signal samples, which limit 
the statistical precision of the form- factor measurement. 



II. FORM FACTORS 



B. Form Factors 



A. Overview 



B Decays to Pseudoscalar Mesons: B — > 1x11/ 



The advantage of charmless semileptonic decays over 
charmless hadronic decays of the B meson is that the 
leptonic and hadronic components of the matrix element 
factorize. The hadronic matrix element is difficult to cal- 
culate, since it must take into account physical mesons, 
rather than free quarks. Therefore higher-order per- 
turbative corrections and non-pcrturbative long-distance 
hadronization processes cannot be ignored. To overcome 
these difficulties, a set of Lorentz-invariant form factors 
has been introduced that give a global description of 
these QCD processes. 

A variety of theoretical predictions for these form fac- 
tors exist. They are based on QCD calculations, such as 
lattice QCD and sum rules, in addition to quark mod- 
els. We will make use of a variety of these calculations 
to assess their impact on the determination of \Vub\ from 
measurements of the decay rates. 

The V — A structure of the hadronic current is invoked, 
along with the knowledge of the transformation proper- 
ties of the final-state meson, to formulate these form fac- 
tors. They arc functions of cf' ~ rn^, the mass squared 
of the virtual W ., 

= iPB-Pxf = Ml+ml-2MBEx. (1) 

Here P( and P^, refer to the four-momenta of the charged 
lepton and the neutrino, Mb and Pb to the mass and 
the four-momentum of the B meson, and mx and Ex 
are the mass and energy (in the B-meson rest frame) of 
the final-state meson X^. 

We distinguish two main categories of exclusive 
semileptonic decays: decays to pseudoscalar mesons, 
B —>■ irii/ or 5+ —>■ ■qi'^iy, and decays to vector mesons, 
B -^ plv or B+ -> oj^+u. 



For decays to a final-state pseudoscalar meson, the 
hadronic matrix clement is usually written in terms of 
two form factors, f+{q^) and fo{q'^) [l3.ll3|. 



{TT{P^)\urb\B{PB)) = 



Ml') 






{Pb + P.Y 

t2 



„9 y 



q- 



(2) 



where P-^ and Pb arc the four-momenta of the final-state 
pion and the parent B meson, and ttItt and AIb arc their 
masses. This expression can be simplified for leptons 
with small masses, such as electrons and muons, because 
in the limit of mi <^ Mb the second term can be ne- 
glected. We are left with a single form factor f+{q^) and 
the differential decay rate becomes 



dr(BO 



-i+v) 



dq^d cos 6 



I Kb 



wi 



327r3 



sm'ewAf+iq'r, (3) 



where p^^ is the momentum of the pion in the rest frame 
of the B meson, and q^ varies from zero to q'^ax = {Mb — 

The decay rate depends on the third power of the pion 
momentum, suppressing the rate at high q^. The rate 
also depends on sin Owii where 9wi is the angle of the 
charged-lepton momentum in the W rest frame with re- 
spect to direction of the W boost from the B rest frame. 
The combination of these two factors leads to a lepton- 
momentum spectrum that is peaked well below the kine- 
matic limit (see Figure [T|). 



2. B Decays to Vector Mesons: B — > plv 
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FIG. 1: Simulated distributions of (f versus En for a) _B — >■ 
■nlv and b) B — > plv decays. Ei is the lepton energy in the 
B-meson rest frame. 

Figure [T] shows the phase space for B -^ irtv and 
B — > pHv decays in terms of q^ and E^, the energy of 
the charged lepton in the i?-meson rest frame. The dif- 
ference between the distributions is due to the different 
spin structure of the decays. 



For decays with a vector meson in the final state, the 
polarization vector e of the vector meson plays an im- 
portant role. The hadronic current is written in terms of 
four form factors, of which only three {Ai with i = 0, 1, 2) 
are independent [12, [3 ' 



pai^BP 



{p{Pp,e)\V^ ^ A^\B{Pb)) = 2»^(g') e^-^e^Pp^P, 

Mb + mp 

-{MB+mp)A^(q^)e*^ + -^^^e* ■ Pb (Pb+PpT 



MB+nip 
+2mp'—^q'^[A,{q') - Ao{q% 



(4) 



where mp and Pp refer to the vcctor-mcson mass and 
four-momentum. Again, a simplification can be made 
for low-mass charged leptons. The term with q^ can be 



neglected, so there are effectively only three form factors 
for electrons and muons: the axial-vector form factors, 
Ai{q'^) and A2{q'^), and the vector form factor, y(g^). 
Instead of using these form factors, the full differential 
decay rate is usually expressed in terms of the helicity 
amplitudes corresponding to the three helicity states of 
the p meson. 



Figure [2] shows the q^ distributions for B — >■ irtiy and 
B — > pif decays for various form-factor calculations. The 
uncertainties in these predictions are not indicated. For 
B — > nil/ decays they are largest at low q^ for LQCD 
predictions and largest at high q^ for LCSR calculations. 
Estimates of the uncertainties of the calculations are cur- 
rently not available for B -^ p(.v decays. 



H±[q^) = [Mb + nip) 



Ai{q')T 



2MbPp 



2rap\f(p- 
{Mb -VvTip)' 



{Mb + rapY 
{Ml-ml-qMii'f) 

:A2{q^) 



V{q^ 



(5) 



where Pp is the momentum of the final-state p meson in 
the B rest frame. While Ai dominates the three helicity 
amplitudes, A2 contributes only to Hq, and V contributes 
only to H±. 

Thus the differential decay rate can be written as 



dT{B ^ piu) ,2 G|Pp9l^ 

-7^- = \V.ub\ -^—^r^^X 

dq^ a COS Owe. 



1287r3M| 



SmOlrAHn^'^ 



7W'£|^0| 



-t-(l — COS^VK^) 



2|ff+ 



-I- (1 + cos6we.) 



2\H-\ 



(6) 



The V—A nature of the charged weak current leads to 
a dominant contribution from _ff_ and a distribution of 
events characterized by a forward peak in cos^vy^ and 
high lepton momenta (see Figure [T]). 



C. Form-Factor Calculations and Models 

The q^ dependence of the form factors can be extracted 
from the data. Since the differential decay rates are pro- 
portional to the product of \Vub\'^ and the form- factor 
terms, we need at least one point in q^ at which the form 
factor is predicted in order to extract \Vub\ from the mea- 
sured branching fractions. 

Currently predictions of form factors are based on 

• quark-model calculations, (ISGW2) [ij], 

• QCD light-cone sum rules (LCSR) [IMll, 



• lattice QCD calculations (LQCD) 120|-l23[. 

These calculations will also be used to simulate the kine- 
matics of the signal decay modes and thus might im- 
pact the detection efficiency and thereby the branching- 
fraction measurement. The two QCD calculations result 
in predictions for different regions of phase space. The 
lattice calculations are only available in the high-g^ re- 
gion, while LCSR provide information near q^ = 0. Inter- 
polations between these two regions can be constrained 
by unitarity and analyticity requirements [2J, |25| . 
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FIG. 2: q^ distributions for B -s> -kIv (left) and B -)■ piv 
(right) decays, based on form-factor predictions from the 
ISGW2 model [Hi, LCSR calculations (LCSR 1 [H and 
LCSR 2 [11 for B ^ ixlv and LCSR [i3| for B ^ plv) 
and the HPQCD [231 lattice calculation. The extrapolations 
of the QCD predictions to the full q range are marked as 
dashed lines. 

The Isgur-Scora-Grinstcin-Wise model (ISGW2) [H 
is a constituent quark model with relativistic corrections. 
Predictions extend over the full q^ range; they are nor- 



malized at q^ 
ized as 



The form factors are parameter- 



/+(g^) = /(94,J(l + — e^(g: 



~imax 



-N 



(7) 



where ^ is the charge radius of the final-state meson, 
and A^ = 2 (A^ = 3) for decays to pseudoscalar (vec- 
tor) mesons. The uncertainties of the predictions by this 
model are difficult to quantify. 

QCD light-cone sum-rule calculations are non- 
perturbative and combine the idea of QCD sum rules 
with twist expansions performed to 0{as). These calcu- 
lations provide estimates of various form factors at low 
to intermediate g^, for both pseudoscalar and vector de- 
cays. The overall normalization is predicted at low q^ 
with typical uncertainties of 10-13% p^[l7|. 

Lattice QCD calculations can potentially provide 
heavy-to-light-quark form factors from first principles. 
Unquenched lattice calculations, in which quark-loop ef- 
fects in the QCD vacuum are incorporated, are now 
available for the B — > ttIv form factors from the Fer- 
milab/MILC [H] and the HPQCD [H] Collaborations. 
Both calculations account for three dynamical quark fla- 
vors, the mass-degenerate u and d quarks and a heavier 
s quark, but they differ in the way the h quark is sim- 
ulated. Predictions for fo{q^) and f+{q^) are shown in 
Figure [3] The two lattice calculations agree within the 
stated uncertainties, which are significantly smaller than 
those of earlier quenched approximations. 
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FIG. 3: (color online) Predictions of the form factors /+ and 
/o for B — >■ nil/ decays based on unquenched LQCD calcula- 
tions by the FNAL/MILC ^ and HPQCD ^ Collabora- 
tions (data points with combined statistical and systematic 
errors) and LCSR calculations [3| (solid black lines). The 
dashed lines indicate the extrapolations of the LCSR predic- 
tions to q^ > 16GeV^. 



D. Form-Factor Parameterizations 

Neither the lattice nor the LCSR QCD calculations 
predict the form factors over the full q^ range. Lattice 
calculations are restricted to small hadron momenta, i.e., 
to q"^ > g^Q^/2, while LCSR work best at small q'^. If 
the q'^ spectrum is well measured, the shape of the form 
factors can be constrained, and the QCD calculations 
provide the normalization necessary to determine \Vub\- 

A number of parameterizations of the pseudoscalar 
form factor f+{q^) are available in the literature. The 
following four parameterizations are commonly used. All 
of them include at least one pole term at q^ ~ "^s* ■, with 
ruB* = 5.325 GeV < Mb + m^. 



1. Becirevic-Kaidalov (BK) [2e 






/+(0) 



(1 - q'^/m%,){l - aBKq'^/m%.)' 
/o(0) 



1 - PsW/rnl, 



(8) 
(9) 



where /+(0) and /o(0) set the normalizations and 
aBK and Pbk define the shapes. The BK param- 
eterization has been applied in fits to the HPQCD 
lattice predictions for form factors, with the con- 
straint/+(0) = /o(0). 

2. Ball-Zwicky (BZ) [H, [T^ : 



/+(?') = ./+(0) 



1 



^l, 



- q'^/rnj 
rBzq^/m^B, 



(1 - q'^/m%,){l - aBzq^/m%,) 



(10) 



where /-i-(0) is the normalization, and aBZ and 
fBZ determine the shape. This is an extension 
of the BK ansatz, related by the simplification 
aBK ~ oiBZ = ^BZ- This ansatz was used to ex- 
tend the LCSR predictions to higher q^, as shown 
in Figure [H 

3. Boyd, Grinstein, Lebed (BGL) ^ [H : 



U{q') 



1 



V{q^)^{q',ql) ^^ 



Y,au{ql)[z{q\ql)]\ (11) 



^(r,9o) = 



ql 



q^ 
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where m± = Mb ± rriT^ and q^ is a free param- 
eter [13. The so-called Blaschke factor V{q^) = 
z{q^,w?g,) accounts for the pole at q^ = m'j,, 
and (j){q^,qQ) is an arbitrary analytic function [23| 
whose choice only affects the particular values of 
the series coefficients Ofc. In this expansion in the 
variable z, the shape is given by the values of Ofc, 
with truncation at k„iax = 2 or 3. The expansion 
parameters arc constrained by unitarity, X]fe '^fc — 
1. Bccher and Hill [23 have pointed out that due 
to the large &-quark mass, this bound is far from 
being saturated. Assuming that the ratio A/mi, is 
less than 0.1, the heavy-quark bound is approxi- 
mately 30 times more constraining than the bound 
from unitarity alone, J2k°'k ~ (^/wfa)^ ~ 0.001. 
For more details we refer to the literature 123, ^ . 



4. Bourrely, Caprini, Lellouch (BCL) [2£ 
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where the variable z is defined as in Eq.[T^with free 
parameter gp (27|. In this expansion the shape is 
given by the values of bk, with truncation at kmax = 
2 or 3. The BCL parameterization exhibits the 
QCD scaling behavior f+{q^) oc 1/q^ at large q^. 

The BK and BZ parameterizations arc intuitive and 
have few free parameters. Fits to the previous BaBAR 
form-factor measurements using these parameterizations 
have shown that they describe the data quite well [^. 
The BGL and BCL parameterizations are based on fun- 
damental theoretical concepts like analyticity and unitar- 
ity. The z-expansion avoids ad hoc assumptions about 
the number of poles and pole masses, and it can be 
adapted to the precision of the data. 



III. DATA SAMPLE, DETECTOR, AND 
SIMULATION 

A. Data Sample 

The data used in this analysis were recorded with the 
BaBAR detector at the PEP-II energy-asymmetric e~^e~ 
coUider operating at the T(4S') resonance. A sample of 
377 million T{AS) — ?► BB events, corresponding to an 
integrated luminosity of 349 fb~^, was collected. An ad- 
ditional sample of 35.1 fb~^ was recorded at a center- 
of-mass (cm.) energy approximately 40 MeV below the 
T(45') resonance, i.e., just below the threshold for BB 
production. This off-resonance data sample is used to 
subtract the non-BB contributions from the data col- 
lected at the T{4:S) resonance. The principal source of 
these hadronic non-BB events is e+e^ annihilation in 
the continuum to qq pairs, where q = u,d,s,c refers to 
quarks. The relative normalization of the off-resonance 
and on-resonance data samples is derived from lumi- 
nosity measurements, which are based on the number 
of detected jjl^ ii~ pairs and the QED cross section for 
e+e^ -^ fj.^fi'^ production, adjusted for the small differ- 
ence in cm. energy. The systematic error on the relative 
normalization is estimated to be 0.25%. 



B. BABAR Detector 

The BaBAR detector and event reconstruction are de- 
scribed in detail elsewhere [30, Isil- The momenta and 
angles of charged particles are measured in a tracking sys- 
tem consisting of a five-layer silicon vertex tracker (SVT) 
and a 40-layer drift chamber (DCH) filled with a helium- 
isobutane gas mixture. Charged particles of different 
masses are distinguished by their ionization energy loss in 
the tracking devices and by a ring-imaging Cerenkov de- 
tector (DIRC). Electromagnetic showers from electrons 
and photons are measured in a finely segmented CsI(Tl) 
calorimeter (EMC). These detector components are em- 
bedded in the 1.5-T magnetic field of the solenoid. The 
magnet flux return steel is segmented and instrumented 
(IFR) with planar resistive plate chambers and limited 
streamer tubes, which detect particles penetrating the 
magnet coil and steel. 

The efficiency for the reconstruction of charged par- 
ticles inside the fiducial volume of the tracking system 
exceeds 96% and is well reproduced by MC simula- 
tion. An effort has been made to minimize fake charged 
tracks, caused by multiple counting of a single low-energy 
track curling in the DCH, split tracks, or background- 
generated tracks. The average uncertainty in the track- 
reconstruction efficiency is estimated to range from 0.25% 
to 0.5% per track. 

To remove beam generated background and noise in 
the EMC, photon candidates are required to have an en- 
ergy of more than 50 MeV and a shower shape that is 
consistent with an electromagnatic shower. The photon 



efficiency and its uncertainty are evaluated by compar- 
ing T -> TT 1/ to T — >■ pi/ samples and by studying 
e+e~ — > /i+/i^(7) events. 

Electron candidates are selected on the basis of the 
ratio of the energy detected in the EMC and the track 
momentum, the EMC shower shape, the energy loss in 
the SVT and DCH, and the angle of the Cerenkov pho- 
tons reconstructed in the DIRC. The energy of electrons 
is corrected for bremsstrahlung detected as photons emit- 
ted close to the electron direction. Muons are identified 
by using a neural network that combines the information 
from the IFR with the measured track momentum and 
the energy deposition in the EMC. 

The electron and muon identification efficiencies and 
the probabilities to misidentify a pion, kaon, or proton as 
an electron or muon are measured as a function of the lab- 
oratory momentum and angles using high-purity samples 
of particles selected from data. These measurements are 
performed separately for positive and negative leptons. 
For the determination of misidentification probabilities, 
knowledge of the inclusive momentum spectra of positive 
and negative hadrons, and the measured fractions of pi- 
ous, kaons and protons and their misidentification rates 
is used. 

Within the acceptance of the SVT, DCH and EMC de- 
fined by the polar angle in the laboratory frame, —0.72 < 
cos^iab < 0.92, the average electron efficiency for labora- 
tory momenta above 0.5 GeV is 93%, largely independent 
of momentum. The average hadron misidentification rate 
is less than 0.2%. Within the same polar-angle accep- 
tance, the average muon efficiency rises with laboratory 
momentum to reach a plateau of about 70% above 1.4 
GeV. The muon efficiency varies between 50% and 80% 
as a function of the polar angle. The average hadron 
misidentification rate is 2.5%, varying by about 1% as a 
function of momentum and polar angle. 

Neutral pious are reconstructed from pairs of photon 
candidates that are detected in the EMC and assumed to 
originate from the primary vertex. Photon pairs with an 
invariant mass within 17.5 MeV of the nominal tt*' mass 
are considered tt*^ candidates. The overall detection effi- 
ciency, including solid angle restrictions, varies between 
55% and 65% for 7r° energies in the range of 0.2 to 2.5 
GeV. 



C. Monte Carlo Simulation 

We assume that the T(4S') resonance decays exclu- 
sively to BB pairs [33 and that the non-resonant cross 
section for e+e^ ^- qq is 3.4 nb, compared to the T{AS) 
peak cross section of 1.05 nb. Wc use Monte Carlo (MC) 
techniques [33i] to simulate the production and decay of 
BB and qq pairs and the detector response [3J], to esti- 
mate signal and background efficiencies, and to extract 
the expected signal and background distributions. The 
size of the simulated sample of generic BB events exceeds 
the BB data sample by about a factor of three, while the 



MC samples for inclusive and exclusive B — >■ Xuli^ de- 
cays exceed the data samples by factors of 15 or larger. 
The MC sample for qq events is comparable in size to the 
qq data sample recorded at the T{4S) resonance. 

Information extracted from studies of selected data 
control samples on efficiencies and resolution is used to 
improve the accuracy of the simulation. Specifically, 
comparisons of data with the MC simulations reveal 
small differences in the tracking efficiencies and calorime- 
ter resolution. We apply corrections to account for these 
differences. The MC simulations include radiative ef- 
fects such as bremsstrahlung in the detector material and 
initial-state and final-state radiation [3a]- Adjustments 
are made to take into account the small variations of the 
beam energies over time. 

For this analysis, no attempt is made to reconstruct 
K^ interacting in the EMC or IFR. Since a K^ deposits 
only a small fraction of its energy in the EMC, K^ pro- 
duction can have a significant impact on the energy and 
momentum balance of the whole event and thereby the 
neutrino reconstruction. It is therefore important to ver- 
ify that the production rate of neutral kaons and their 
interactions in the detector are well reproduced. 

From detailed studies of large data and MC samples 
of Z?" — > i^r^TT+TT" and Z?" — >■ T^^tt+tt" decays, correc- 
tions to the simulation of the K^ detection efficiency and 
energy deposition in the EMC are determined. The MC 
simulation reproduces the efficiencies well for Kl labo- 
ratory momenta above 0.7 GeV. At lower momenta, the 
difference between MC and data increases significantly; 
in this range the MC efficiencies are reduced by randomly 
eliminating a fraction of the associated EMC showers. 
The energy deposited by K^ in the EMC is significantly 
underestimated by the simulation for momenta up to 
1.5 GeV. At higher momenta the differences decrease. 
Thus the simulated energies are scaled by factors vary- 
ing between 1.20 and 1.05 as a function of momentum. 
Furthermore, assuming equal inclusive production rates 
for K^l and Kg we verify the production rate as a function 
of momentum, by comparing data and MC simulated Kg 
momentum spectra. We observe differences at small mo- 
menta; below 0.4 GeV the data rate is lower by as much 
as 22 ± 7% compared to the MC simulation. To account 
for this difference, we reduce the rate of low momentum 
K1 in the simulation by randomly transforming the ex- 
cess A'2 into a fake photon, i.e., we replace the energy 
deposited in the EMC by the total K^^ energy and set 
the mass to zero. Thus we correct the overall energy 
imbalance created by the excess in K'^ production. 

For reference, the values of the branching fractions, 
lifetimes, and parameters most relevant to the MC sim- 
ulation are presented in Tables HI and HIl 

The simulation of inclusive charmless semileptonic de- 
cays B — >■ Xuiv is based on predictions of a heavy-quark 
expansion (HQE) (valid to 0{as) jSg) for the differen- 
tial decay rates. This calculation produces a smooth 
hadronic mass spectrum. The hadronization of X„ with 
masses above 2?7i^ is performed by JETSET [33] ■ To 



TABLE I: Branching fractions and their errors for the semilep- 
tonic B decays used in this analysis. 



Decay 


Unit 


B° B± 


Ref. 


B-i> rjeu 


10-* 


0.40 ± 0.09 


m 


B -5> Tj'iu 


10"* 


0.21 ±0.21 


m 


B ^iulu 


10-* 


1.15 ±0.16 


m 


B -^ Xuiv 


10"^ 


2.25 ±0.22 2.41 ±0.22 


m 


B^ Dlu 


10-2 


2.17 ±0.08 2.32 ±0.09 


m 


B^ D*lu 


10-2 


5.11 ±0.19 5.48 ±0.27 


[411 


B -^ Dllv 


10-2 


0.69 ±0.14 0.77 ±0.15 


[41,45] 


B -^ D*2lu 


10-2 


0.56 ±0.11 0.59 ±0.12 


[41,45] 


B -^ DU^ 


10-2 


0.81 ±0.24 0.88 ±0.26 


[41,45] 


B -^ D'Ju 


10-2 


0.76 ±0.22 0.82 ±0.25 


[41,45] 



TABLE II: Form factors used in the simulation of B — >■ Div 
and B — > D* iv decays, based on the parameterization of 
Caprini, Lellouch and Neubert ^^. The B° lifetime, the _B° 
to B^ lifetime ratio, and relative branching fraction at the 
T(4S) resonance, used in the simulation. 



Parameter 


Value 


Rcf. 


B ^ Dlu FF ■ pi, 


1.18 ±0.04 ±0.04 


m 


B-^ D'lvFF ■ pi. 


1.191 ±0.048 ±0.028 


m 


B ^ D'luFF : Ri 


1.429 ± 0.061 ± 0.044 


m 


B^ D'lvFF:R2 


0.827 ±0.038 ±0.022 


m 


B° life time to ( ps) 


1.530 ± 0.009 


m 


Blife time ratio t+/tq 


1.071 ± 0.009 


m 


r(4S) ratio /+_//oo 


1.065 ± 0.026 


[41] 



describe the dynamics of the h quark inside the B me- 
son we use HQE parameters extracted from global fits 
to moments of inclusive Icpton-energy and hadron-mass 
distributions in _B — >■ Xcf-v decays and moments of in- 
clusi ve p hoton-energy distributions in B — >■ Xg^ de- 
cays [S^ . The specific values of the HQE parameters in 
the shape-function scheme are mi, = 4.631 ± 0.034 GeV 
and jjL^ = 0.184 ± 0.36 GeV ; they have a correlation of 
p = —0.27. Samples of exclusive semileptonic decays 
involving low-mass charmless mesons (tt, p, w, ry,??') are 
simulated separately and then combined with samples of 
decays to non-resonant and higher-mass resonant states, 
so that the cumulative distributions of the hadron mass, 
the momentum transfer squared, and the lepton momen- 
tum reproduce the HQE predictions. The generated dis- 
tributions arc reweighted to accommodate variations due 
to specific choices of the parameters for the inclusive and 
exclusive decays. The overall normalization is adjusted 
to reproduce the measured inclusive B — )■ X^iv branch- 
ing fraction. 

For the generation of decays involving charmless 
pseudo-scalar mesons we choose two approaches. For the 
signal decay B — )■ -kIv we use the ansatz by Becirevic 
and Kaidalov [2g for the q^ dependence, with the single 
parameter ubk set to the value determined in a previous 



BaBAR analysis Q of B -^ 7r£i/ decays, asK = 0.52±0.06. 
For decays to rj and rj' we use the form factor parameter- 
ization of Ball and Zwicky with specific values reported 
in [3. 

Decays involving charmless vector mesons (p, 1^) arc 
generated based on form factors determined from LCSR 
by Ball and Zwicky |17| . We use the parameterization 
proposed by the authors to describe the q^ dependence 
of the form factors in terms of a modified pole ansatz us- 
ing up to three independent parameters ri, r2 and mm- 
Table IIIII shows the suggested values for these parame- 
ters. TTifit refers to an effective pole mass that accounts for 
contributions from higher- mass B mesons with J^ ==1", 
and ri, and r2 give the relative scale of the two pole 
terms. 



TABLE III: Parameterization of the LCSR form-factor calcu- 
lations [i^, [131 for decays to pseudo-scalar mesons 77 and rj' 
(/+) and vector mesons p and to {Ai, A2,V). 



Form factor /+ A^ 



A', 



V 



^1 A^ 



1/- 



F(0) 0.273 0.242 0.221 0.323 0.219 0.198 0.293 

n 0.122 - 0.009 1.045 - 0.006 1.006 

r2 0.155 0.240 0.212 -0.721 0.217 0.192 -0.713 

mlt(GeV^) 31.46 37.51 40.82 38.34 37.01 41.24 37.45 
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ing fractions by 60% and inflate the errors by a factor of 
three. 



IV. EVENT RECONSTRUCTION AND 
CANDIDATE SELECTION 

In the following, we describe the selection and kine- 
matic reconstruction of signal candidates, the definition 
of the various background classes, and the application of 
neural networks to suppress these backgrounds. 



A. Signal-Candidate Selection 

Signal candidates are selected from events having at 
least four charged tracks. The reconstruction of the 
four signal decay modes, B^ — > tt~£^i^, i?+ -> 7r°£+z/, 
B'^ -^ p~t^v and B^ — > p^t^v^ requires the identi- 
fication of a charged lepton, the reconstruction of the 
hadronic state consisting of one or more charged or neu- 
tral pions, and the reconstruction of the neutrino from 
the missing energy and missing momentum of the whole 
event. 



Lepton and Hadron Selection 



For the simulation of the dominant B — > Xciv decays, 
we have chosen a variety of models. For B — > Dlv and 
B — >■ 0*11/ decays we use parameterizations [33, EO] of 
the form factors based on heavy quark effective theory 
(HQET). In the limit of negligible lepton masses, decays 
to pseudoscalar mesons are described by a single form 
factor for which the q^ dependence is given by a slope 



parameter. We use the world average 14]|, updated for 
recent precise measurements by the BaBAR Collabora- 
tion [42, I43 . Decays to vector mesons are described by 
three form factors, of which the axial vector form fac- 
tor dominates. In the limit of heavy quark symmetry, 
their q^ dependence can be described by three parame- 
ters: p|)», i?i, and i?2. We use the most precise BABAR 
measurement [44| of these parameters. 

For the generation of the semileptonic decays to D** 
resonances (four L = 1 states), we use the ISGW2 [lj| 
model. At present, the sum of the branching fractions 
for these four decays modes is measured to be 1.7%, but 
so far only the decays D** —)■ Dn and D** — > D*Tr have 
been reconstructed, while the total individual branching 
fractions for these four states remain unknown. Since the 
measured inclusive branching fraction for B — > Xcii' ex- 
ceeds the sum of the measured branching fractions of all 
exclusive semileptonic decays by about 1.0%, and since 
non- resonant B — > D'^*^'kI!.v decays have not been ob- 
served [4^, we assume that the missing decays are due 
to i3 — >■ D**£v, involving hadronic decays of the D** 
mesons that have not yet been measured. To account for 
the observed deficit, we increase the B —J- D**£v branch- 



Candidates for leptons, both e and p^ , are required 
to have high cm. momenta, p*^ > 1.0 GeV for the B — J- 
7t£i^, and p| > 1.8 GcV for the B — > pli/ sample. This 
requirement significantly reduces the background from 
hadrons that are misidentified as leptons, and also re- 
moves a large fraction of true leptons from secondary 
decays or photon conversions, and from B — > Xciv de- 
cays. 

To suppress Bhabha scattering and two-photon pro- 
cesses in which an electron or a photon from initial- 
state or final state radiation interacts in the material of 
the detector and generates additional charged tracks and 
photons at small angles to the beam axis, we require 
^z < 0.65 for events with a candidate electron. Here 
?2 ~ J2iPi/ J2i ^i' where the sum runs over all charged 
particles in the event and pf and Ei are their longitu- 
dinal momentum components and energies measured in 
the laboratory frame. 

For the reconstruction of the signal hadron, we con- 
sider all charged tracks that are not consistent with a 
signal lepton and not identified as a kaon. Neutral pions 
are reconstructed from pairs of photon candidates and 
the 7r° cm. momentum is required to exceed 0.2 GcV. 
Candidate p^ — >■ tt^tt" or p° — > tt+tt^ decays are re- 
quired to have a two-pion mass within one full width of 
the nominal p mass, 0.650 < A/^tt < 0.850 GeV. To re- 
duce the combinatorial background, we also require that 
the cm. momentum of one of the pions exceed 0.4 GeV, 
and that the cm. momentum of the other pion be larger 
than 0.2 GeV. 
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Each charged lepton candidate is combined with a 
hadron candidate to form a so-called Y candidate of 
charge zero or one. At this stage in the analysis we allow 
for more than one candidate per event. Two or three 
charged tracks associated with the Y candidate are fit- 
ted to a common vertex. This vertex fit must yield a 
X^ probability of at least 0.1%. To remove background 
from J/i/; —>■ ^^i,~ decays, we reject a Y candidate if the 
invariant mass of the lepton and any oppositely charged 
track in the event is consistent with this decay. 

To further reduce backgrounds without significant sig- 
nal losses, we impose additional restrictions on the cm. 
momenta of the lepton and hadron candidates by requir- 
ing at least one of the following conditions to be satisfied, 
for B — > tt£v 





pL, 


,,„„>1.3GeV 


or 






p*f > 2.2 GcV 


or 




PLd™„+rf>2.8GeV, 


iOT B 


-^ pHv 








pL 


,,,„„> 1.3 GeV 


or 






rf>2.0GeV 


or 




pLdron+P*i>2.65GeV. 



These additional requirements on the lepton and hadron 
cm. momenta primarily reject background candidates 
that are inconsistent with the phase space of the signal 
decay modes. 

If a K candidate originates from a signal decay mode, 
the cosine of the angle between the momentum vectors 
of the B meson and the Y candidate, cos 6'By, can be 
calculated as follows. 



cos t'BY 



zEj^Ey 



Ml 



M2 



2p*g pl- 



(14) 



and the condition |cos(?by| < 1-0 should be fulfilled. 
The energy Eg and momentum p|j of the B meson are 
not measured event-by-event. Specifically, Eg = \/s/2 is 
given by the average cm. energy of the colliding beams, 
and the B momentum is derived as p*g = \/ E*b ^ ^^b- 
To allow for the finite resolution in this variable, we im- 
pose the requirement —1.2 < cos6'by < l-l- 



2. Neutrino Reconstruction 



frame. The energy calculation depends on the correct 
mass assignments for charged tracks. For this reason we 
choose to calculate the missing momentum and energy in 
the laboratory frame rather than in the rest frame of the 
TiAS). By doing so, we keep this uncertainty confined 
to the missing energy. 

If all particles in the event, except the neutrino, are 
well measured, P^ ~ (£'miss,Pmiss) is a good approxima- 
tion. However, particles that are undetected because of 
inefficiency or acceptance losses, in particular Ki^ mesons 
and additional neutrinos, or spurious tracks or photons 
that do not originate from the BB event, impact the 
accuracy of this approximation. To reduce the effect of 
losses due to the limited detector acceptance, we require 
that the polar angle of the missing momentum in the 
laboratory frame be in the range 0.3 < ^miss < 2.2 rad. 
We also require the missing momentum in the laboratory 
frame to exceed 0.5 GeV. 

For the rejection of background events and signal de- 
cays that are poorly reconstructed as well as events with 
more than one missing particle, we make use of the miss- 
ing mass squared of the whole event. 



P^ 



Ei 



\Pn 



(16) 



For a correctly reconstructed event with a single semilep- 
tonic B decay, JTi^jsg should be consistent with zero 
within measurement errors. Failure to detect one 
or more particles in the event creates a substantial 
tail at large positive values. Since the resolution in 
TTi^jgg increases linearly with i^miss, we use the variable 
TO^igg/2i?iniss — E^iss - Pmiss as a discriminator and re- 
quire m'^^j2En^iss < 2.5 GeV. 



3. Variables Used for Signal Extraction 

The kinematic consistency of the candidate decay with 
a signal B decay is ascertained using two variables, the 
beam-energy substituted B mass toes i and the difference 
between the reconstructed and expected energy of the B 
candidate AE. In the laboratory frame, they are defined 
as 



AE- 



Pb ■ Pe 



s/2 



^s 



(17) 



and 



The neutrino four- momentum, P^, ~ (^miss, Pmiss), 
is inferred from the difference between the net four- 
momentum of the coUiding-beam particles, Pe+e- = 
(i?e+e- jPe+e- ), and the sum of the measured four-vectors 
of all detected particles in the event, 

(Smiss, Pmiss) = (^e+e" , Pe+e" ) - C^^i.^p^), (15) 



where Ei and pi are the energy and three-momentum of 
the i"^ track or EMC shower, measured in the laboratory 



TOES 



'{s/2+pB -Pe+e-)' 



E' 



Pi 



(18) 



where Pb = {Eb,Pb) and Pg+e- denote the four- 
momenta of the B meson and the colliding beam par- 
ticles, respectively. The B-meson momentum vector pB 
is determined from the measured three-momenta of the 
decay products, and Pg+e- is derived from the calibra- 
tion and monitoring of the energies and angles of the 
stored beams. We extract the signal yields by a fit to the 
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two-dimensional AE — tties distributions in bins of the 
momentum transfer squared q^. We define a region in 
the AE — rriEs plane that contains almost all of the sig- 
nal events and leaves sufficient phase space to constrain 
the different background contributions. This fit region is 
defined as 

I A£;| < 0.95 GeV , 5.095 < toes < 5.295 GeV. (19) 

Only candidates that fall inside the fit region are con- 
sidered in the analysis. We also define a smaller region 
where the signal contribution is much enhanced relative 
to the background. This signal region is defined as 

- 0.15 <AE< 0.25 GeV , 5.255 < ?«es < 5.295 GeV. 

.(20) 
The signal region is chosen to be slightly asymmetric in 
AE to avoid sizable B — > Xcii' background, which peaks 
near —0.2 GeV. In the following, we refer to the phase 
space outside the signal region, but inside the fit region, 
as the side bands. 

As a measure of the momentum transfer squared q^ 
we adopt the mass squared of the virtual W, i.e., the 
invariant mass squared of the four-vector sum of the re- 
constructed lepton and neutrino. 



lu 



[iEi,Pi) + {Eyniss,Pmiss)Y 



(21) 



The resolution in q^^.^ is dominated by the measurement 
of the missing energy which tends to have a poorer res- 
olution than the measured missing momentum, because 
the missing momentum is a vector sum and contributions 
from particle losses (or additional tracks and EMC show- 
ers) do not add linearly as is the case for i^miss- Thus for 
the definition of q'^^^ it is advantageous to replace i^miss 
by Pmiss, the absolute value of the measured missing mo- 
mentum, 



^raw 



[(Ei.pi) + {Pmiss,Pmiss)Y 



(22) 



The resolution of q'^^^ can be further improved by scal- 
ing Pmiss by a factor of a, such that AE of the B candi- 
date is forced to zero, 



p^ = apn 



with 



1 



AE 



(23) 



and substituting pi^ for Pmiss to obtain gcon- ^^ly candi- 
dates for which this q^ correction yields unphysical val- 
ues, g^Qj,^ < OGeV^, are rejected. This is the case for 
about 1% of the background not associated with semilep- 
tonic decays. The quantity (/^orr is used as the measured 
q^ throughout this analysis. 

The q^ resolution is critical for the measurements of the 
form factors. Figure 2] shows the correlation between the 
true q^ and the reconstructed g^o^r fo^' simulated sam- 
ples of B — )- tt£i' and B — >■ piv candidates passing the 
entire event selection, which is described below. Cor- 
rectly reconstructed signal events and combinatorial sig- 
nal events, for which the hadron has been incorrectly 



selected, are shown. For correctly reconstructed signal 
decays the resolution improves with higher q^ and can 
be well described by the sum of two Gaussian resolution 
functions, see Table ITVl In the signal region, the widths 
of the core resolution are in the range 0.18 — 0.34 GeV^, 
and the tails can be approximated by a second Gaus- 
sian function with widths in the range 0.6 — 0.8 GeV^. 
As expected, the resolution is significantly worse in the 
larger fit region. Combinatorial signal events contribute 
primarily at high q^. We rely on the Monte Carlo sim- 
ulation to reproduce the resolution in the reconstructed 
g^Qjj. variable. 



B. Background Suppression 

1. Signal and Background Sources 

A variety of processes contribute to the four samples of 
selected candidates for the charmless semileptonic decay 
modes B° -^ tt-I+i^, B+ -^ n°t+v, B° -^ p't+v, and 
5+ — >■ p^t^v. We divide the signal and background for 
each of the four candidate samples into a set of sources 
based on the origin of the charged lepton candidate. 

• Signal: We differentiate four classes of signal 
events; for all of them the lepton originates from 
a signal decay under study: 

— True signal: the hadron originates from the 
signal decay under study; 

— Combinatorial Signal: the hadron is in- 
correctly selected, in many cases from decay 
products of the second B meson in the event; 

— Isospin-conjugate signal: the lepton origi- 
nates from the isospin conjugate of the signal 
decay; 

— Cross-feed signal: the lepton originates from 
another signal decay mode, for instance B — > 
plv in a _B — >■ T:£iy sample. 

• Continuum background: We differentiate two 
classes of continuum backgrounds: 

— True leptons: the lepton candidate originates 
from a leptonic or semileptonic decay of a 
hadron produced in e+e^ — > qq (mostly cc) 

— > ^"'"£"(7) processes, where £+^~ 



,p+p- 



stands for e'^e 
77; 

— Fake leptons: the lepton candidate is a 
misidentified hadron; this is a sizable contri- 
bution to the muon sample. 

• B ^ Xuiv background: We differentiate two dif- 
ferent sources oi B ^f Xuiv background: 

— Exclusive B — > Xuiv decays involving a sin- 
gle hadron with mass below 1 GeV.' decays 
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FIG. 4: q^ resolution for B" — > -k^I^v (left) and B" — > p^t^v (right) samples after the full event selection: two-dimensional 
distribution of g^ruo versus q^oii in the fit region and in the signal region. Top row: true signal decays, middle row: all 
signal decays (true and combinatorial), bottom row: distribution of g^orr ~ 'Jtruo for true signal (black, solid histogram) and 
combinatorial signal (blue, dashed histogram) in the fit region. The fit of the sum of two Gaussian functions to the true signal 
distribution is shown as a solid red line, the contribution of the broader of the two functions is shown as a dotted red line. 



that are not analyzed as signal (5+ -^ ujl'^iy, 
B+ -^ ryf+i/, and'B+ -^ rfl+i^) ; 

— Inclusive B — > Xuiiy decays: decays involv- 
ing more than one hadron or a single hadron 
with mass above 1 GeV. 

• BB background: We differentiate three classes 
of BB background, excluding B -^ Xuiv decays: 

— Primary leptons, i.e., B — > Xc^v decays: the 
lepton originates from a charm semileptonic 
B decay, either B —>■ Dlv, B — >• D*£h', or 
B -^ D^*^mT)iiy with n> I additional pions; 
this class is dominated by B — >■ D*£i/ decays; 
the largest contributions involve hadrons that 
do not originate from the semileptonic decay; 

— Secondary leptons: the lepton originates 
from the decay of a particle other than a B 



meson, for instance charm mesons, r leptons, 
J/ip , or from photon conversions; 
— Fake leptons: the lepton candidate is not 
a lepton, but a misidcntificd charged hadron; 
this background is dominated by fake niuons. 

Given that the secondary-lcpton and fake-lepton 
BB background contributions are relatively small 
in this analysis, we combine them into one class 
(other BB). 

For intermediate values of q^ (in the range 4 < q^ < 
20GeV^), B — >■ Xcii' decays are by far the dominant 
background, whereas continuum background contributes 
mostly at low and high q^. The B — ?> X„^z/ decays have 
much smaller branching fractions, but their properties 
are very similar to the signal decays and thus they are 
difficult to discriminate against. They contribute mostly 
at high q^, where they are the dominant background. 
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the A_E — niES fit region and in the signal region, integrated over q ; ai, fii and (72, fi2 denote the means and the widths of the 
two Gaussian functions, and the last column lists the fraction of the events characterized by the narrower resolution function. 







Gaussian Fct. 1 


Gaussian Fct. 2 






Signal mode 


Mi(GeV2) 


ai(GeV2) 


/i2(GeV") 


0-2 (GeV^) 


Fraction 


Fit region 


B" ^ TT-i+l. 


-0.005 


0.380 


-0.021 


1.270 


0.35 




B+ ^ TT^'t+U 


0.076 


0.468 


-0.039 


1.343 


0.43 




B» ^p~t+v 


0.005 


0.343 


-0.386 


1.094 


0.45 




B+ ^ p^'l+v 


-0.032 


0.311 


-0.498 


1.086 


0.46 


Signal region 


B" -^ -R-t+u 


0.006 


0.242 


-0.020 


0.720 


0.45 




B+ ^ -kH-^u 


0.058 


0.338 


0.172 


0.807 


0.58 




B° -> p-e.+v 


0.042 


0.246 


0.036 


0.647 


0.50 




B+ ^ pH+u 


0.010 


0.177 


-0.078 


0.586 


0.46 



2. Neural Networks 

To separate signal events from the background sources, 
continuum events, non-signal B — > Xu(.v decays and the 
remaining BB events, we employ a neural-network tech- 
nique based on a multi-layer perceptron (MLP) [46[. We 
have set up a network structure with seven input neu- 
rons and one hidden layer with three neurons and have 
adopted the method introduced by Broyden, Fletcher, 
Goldfarb, and Shanno [43] to train the network. Some 
of the input variables are used as part of the event pre- 
selection that is designed to reduce the BB and contin- 
uum backgrounds by cutting out regions where the signal 
contribution is small or where there are spikes in distri- 
butions, which the neural network may not deal with ef- 
fectively. The following variables are input to the neural 
networks: 

• i?2, the second normalized Fox- Wolfram mo- 
ment |48| determined from all charged and neutral 
particles in the event; we require R2 < 0.5; 

• L2 = X]iPi<^os^6',*, where the sum runs over all 
tracks in the event excluding the Y candidate, and 
p* and 9* refer to the cm. momenta and the angles 
measured with respect to the thrust axis of the Y 
candidate; we set a loose restriction, L2 < 3.0 GeV. 

• cos A6'thrust , where A6'thrust is the angle between 
the thrust axis of the Y candidate and thrust axis 
of all other detected particles in the event; there is 
no preselection requirement for this variable; 



• "i^i,,/(2£;n,iss) = {E'i-^^^ - pf^.J/{2En,i,,); we re- 



l;2 . _ p2 

iniss i-^m: 

< 2.5 GeV; 



quire \mf^^-,^j2E„ 

• cos6'by; we require —1.2 < cos6'by < 1-1; 

• cos 9w£, the helicity angle of the lepton; we require 
|cos6iwd < 0.8 for the B° -^ tt'I+v and B+ ^■ 
'k'^I^v modes; 

• ^missj the polar angle of the missing momentum in 
the laboratory frame; we require 0.3 < O-^-^^ss < 2.2 
rad. 



The first three input variables are sensitive to the topo- 
logical difference between the jet-like continuum events 
and the more spherical BB events. Restrictions on these 
variables do not bias the q^ distribution significantly. 

The restrictions placed on cos^by, "^miss/(2^miss), 
and 0miss do not significantly bias the q^ distribution ei- 
ther. However, the variable cos6'vy^ is correlated with the 
lepton momentum and thereby q^. To ensure that the se- 
lection does not adversely affect the measurement of the 
q^ spectrum, wc have chosen rather moderate restrictions 
on cos 9wi. 

Figure [5] shows the AE and toes distributions for sam- 
ples of i?" — ?► 7r~£+:/ and B'^ — >■ p~i.'^v candidates (inte- 
grated over q^) that have been preselected by the criteria 
described above. The stacked histograms show the sig- 
nal and background contributions compared to the data, 
prior to the fit. The three dominant backgrounds are 
B ^- Xciv decays (including B — >• Div, B — >• D*^^ and 
B -^ D^*'{m:)(.v), qq continuum and B -^ Xuiv decays. 
The signal contributions are very small by comparison 
and difficult to observe. 

The neural networks are trained separately for the 
three background categories and for different q^ intervals. 
Wc introduce six bins in q^ for B — > i:lv and three bins 
for B — > ptu. The bin sizes are 4 GeV^ for B -^ -kIv and 
8 GeV^ for B — )■ plv, except for the last bin, which ex- 
tends to the kinematic limit of 26.4 GeV^ and 20.3 GeV^, 
respectively. Thus in total we train 3 x (2 x 6-1-2 x 3) = 54 
neural networks. Since we aim for a good signal-to- 
background ratio in the region where most of the signal is 
located, we do not train the neural network with events 
in the whole fit region, but in an extended signal region, 
-0.25 < AS < 0.35 GeV, 5.240 < toes < 5.295 GeV. 

For the training of the neural networks we use MC sim- 
ulated events containing correctly reconstructed signal 
decays and the following simulated background samples: 



1. a sample of continuum events, e+e 
gf = u, d, s, c [iqq neural network); 



> qq with 
XcHv decays (B — > 
3. a sample of inclusive B — > X^^-v decays (B — )• X^iv 



2. a combined sample of i? - 
Xciv neural network); and 
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FIG. 5: (color online) Distributions of AE and tties for 
B° -^ ■K^e.+v (left) and _B° -i> p't+v (right) candidates after 
the preselection, i.e., prior to the neural- network application. 
The stacked histograms show the predicted signal and back- 
ground contributions prior to the fit. The expected signal 
distribution (with arbitrary normalization) is indicated as a 
magenta dashed histogram. 



neural network). 

The training of the neural networks and the subsequent 
background reduction is performed sequentially for the 
three background samples. We use subsamples of typi- 
cally less than half the total MC samples for training and 
validation of the neural networks. Of these subsamples, 
one half of the events is used as training sample, and the 
other half for validation. 

Studies of the neural-network performance for the 
B -^ Xu(-v background indicate that the separation of 
this background from the signal is very difficult because 
of the similarity in the shape of the distributions, espe- 
cially for the B'^ — >■ tt^^+z/ and the B -^ plv samples. 
Given these difficulties, we use the B — > Xu^v neural 
network only for the B^ — > 7r~£+i/ sample, and only for 
q^ > 12 GeV , where the B -^ Xu^v background becomes 
significant. 

Figure [6] shows, for the sample of i?° -^ ti~1^v can- 
didates, the distributions of the seven input variables to 
the neural networks. The distributions are shown sequen- 
tially after application of the preselection, the qq neural 
network and the B — >■ Xc^-v neural network to illustrate 
the change in the sample composition. Figures [7] to [S] 
show the distributions of the three neural-network dis- 
criminators for the B'^ — > 'k~ (.^v sample in four of the six 



q^ bins. Figures [TOl and [Til show the distributions of the 
two neural- network discriminators for the B^ — > p~ i~^v 
sample in all three q^ bins. The discriminator cuts are 
chosen to minimize the total error on the signal yield for 
each channel, using the sum in quadrature of the error 
obtained from the maximum-likelihood fit described in 
Section rvTl and the estimated total systematic error of the 
partial signal branching fraction in each q^ bin (see Sec- 
tion |VlTl)- The data-MC agreement is reasonably good 
for the input distributions and the neural-network dis- 
criminators. One should keep in mind that at this stage 
the distributions are taken directly from the simulation, 
without any adjustments or fit. 



TABLE V: Overview of the selection efflciencies for the four 
signal decays (true and combinatorial signal combined) and 
their primary background sources, B — >■ Xulv, B — >■ Xclv, 
and non-BB background. 



Selection 
Units 



Signal Xutv Xdv 

^Q-2 ^Q-3 ;^Q-4 



99 
10-5 



B" 



Preselection 


18.1 


25.6 


26.4 


19.4 


NNgg 


11.6 


13.4 


18.3 


2.6 


NN B -^ Xdv 


7.8 


8.1 


3.4 


1.3 


NN B -^ Xulv 


6.8 


5.3 


2.5 


1.0 


Signal region 


1.8 


0.5 


0.1 


<0.1 


B+ ^ -K^'t^u 










Preselection 


12.8 


20.0 


17.4 


15.1 


NN qq 


8.4 


11.3 


13.3 


2.7 


NN B -^ Xclv 


5.9 


6.6 


1.9 


1.6 


Signal region 


1.6 


0.5 


<0.1 


<0.1 


B" ^ p-l+u 










Preselection 


8.9 


23.9 


35.8 


13.2 
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Table fVl shows the selection efficiencies for the four sig- 
nal samples compared to the efficiencies for the dominant 
background sources for these samples. The total signal 
efficiencies are typically 6 — 7% for B — > niv decays and 
roughly 1 — 2.5% for B — > p(,v decays in the fit region. 
The dominant BB and qq backgrounds are suppressed 
by factors of order 10'' and 10^, respectively. 



3. Candidate Multiplicity 

After the neural-network selection there are on aver- 
age 1.14 candidates per event in the B^ -^ n~i^v sample, 
1.46 in the B+ -^ n°e+v sample, 1.30 in the B° -^ p-£+v 
sample, and 1.17 in the _£?"*" —J- p^t^v sample. We observe 
fewer candidates for decay modes without neutral pious 
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FIG. 6: (color online) Background suppression for B'^ — >■ tt~£'^u candidates. Distributions of the seven input variables to the 
neural network: after the preselection (left column), after the qq neural network (center column), and after the B — >■ Xciv 
neural network (right column). The data are compared to the sum of the MC-simulated background contributions; for a legend 
see Figure [5] The expected signal distribution is overlaid as a magenta, dashed histogram with arbitrary normalization. 
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FIG. 7: (color online) The qq neural-network discriminators for B" — ;> vr i'^iy candidates in the signal region, —0.15 < |A_E| < 
0.25 GeV; 5.255 < tties < 5.295 GeV. The distributions are shown for four different q^ bins, columns from left to right: 
< g^ < 4GeV^, 4 < g^ < 8GeV^ 12 < g^ < 16GeV^ g^ > 20GeV^. Top row: Discriminator distributions for signal 
(magneta, dashed) and gg background (blue, solid), normalized to the same area. The arrows indicate the chosen cuts. Bottom 
row: Discriminator distributions for data compared with MC-simulated signal and background contributions. For a legend see 
Figure [5] 
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5.255 < 



mES 

2 



Xctv neural-network discriminators for B" 
< 



t^v candidates in the signal region, —0.15 < 



5.295 GeV. 



The distributions are shown for four different 
2 ^ icr<„Ar2 2 -. r)nr<„^A2 



bins, columns from left to 



FIG. 

lA^l < 0.25 GeV: 

right: < g^ < 4 GeV, 4 < g^ < 8GeV^, 12 < g^ < 16GeV^, g^ > 20GeV^ Top row: Discriminator distributions for signal 
(magneta, dashed) and B — >■ Xclv background (blue, solid), normalized to the same area. The arrows indicate the chosen cuts. 
Bottom row: Discriminator distributions for data compared with MC-simulated signal and background contributions. For a 
legend see Figure (5] 
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FIG. 9: (color online) The B — )■ Xu£i^ neural-network discriminators for B" — >■ n^l'^u candidates in the signal region, —0.15 < 
|A_E| < 0.25 GeV ; 5.255 < rriES < 5.295 GeV. The distributions are shown for the three highest q^ bins, columns from left 
to right: 12 < g^ < 16 GeV '^j 16 < q^ < 20GeV^, q^ > 20GeV^. Top row: Discriminator distributions for signal (magneta, 
dashed) and B — >■ Xulv background (blue, solid), normalized to the same area. The arrows indicate the chosen cuts. Bottom 
row: Discriminator distributions for data compared with MC-simulated signal and background contributions. For a legend see 
Figure [5] 
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FIG. 10: (color online) The qq neural-network discriminators for 5° — >■ p~ t^v candidates in the signal region, —0.15 < | A^E] < 
0.25 GeV; 5.255 < rriES < 5.295 GeV. The distributions are shown for three different q^ bins, columns from left to right: 
< g'^ < SGeV'^, 8 < g^ < 16 GeV ■^, q^ > 16GeV^. Top row: Discriminator distributions for signal (magneta, dashed) and 
qq background (blue, solid), normalized to the same area. The arrows indicate the chosen cuts. Bottom row: Discriminator 
distributions for data compared with MC-simulated signal and background contributions. For a legend see Figure [5] 



in the final state. For all four samples, the observed can- 
didate multiplicity is well reproduced by MC simulation. 

In case of multiple candidates for a given decay mode, 
we select the one with the highest probability of the ver- 
tex fit for the Y candidate. Since this is not an option for 
i?+ — >■ tt'^I^i' decays, we select the photon pair with an 
invariant mass closest to the tt" mass. Simulations of sig- 
nal events indicate that this procedure selects the correct 
signal decay in 60 — 65% of the cases. By this selection 
we do not allow a single event to contribute more than 
one candidate to a given decay-mode sample, though we 
do allow an event to contribute candidates to more than 



one decay-mode sample. 



V. DATA-MONTE CARLO COMPARISONS 

The determination of the number of signal events relies 
heavily on the MC simulation to correctly describe the 
distributions for signal and background sources. There- 
fore a significant effort has been devoted to detailed com- 
parisons of data and MC distributions for samples that 
have been selected to enhance a given source of back- 
ground. 
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FIG. 11: (color online) The B — >■ Xc£u neural-network discriminators for B" — )> p^t^v candidates in the signal region, 
—0.15 < |A_E| < 0.25 GeV ; 5.255 < mES < 5.295 GeV. The distributions are shown for three different q^ bins, columns from 
left to right: < g^ < 8GeV^, 8 < q^ < 16GeV^, (f > 16GeV^. Top row: Discriminator distributions for signal (magneta, 
dashed) and B — >■ Xclv background (blue, solid), normalized to the same area. The arrows indicate the chosen cuts. Bottom 
row: Discriminator distributions for data compared with MC-simulated signal and background contributions. For a legend see 
Figure [5] 



A. Comparison of Off-Resonance Data with qq MC 

Though we record data below BB threshold (ofF- 
resonance data) the total luminosity of this sample is 
only about 10% of the T{AS) data sample (on-resonance 
data), and thus we need to rely on MC simulation to 
predict the shapes of these background distributions. 

To study the simulation of qq events, we scale the MC 
sample to match the integrated luminosity of the off- 
resonance data. The study is performed separately for 
samples with electrons and muons. This background con- 
tains events with true leptons from leptonic or semilep- 
tonic decays of hadrons, as well as hadrons misidentified 
as leptons. The muon sample is dominated by misidenti- 
fied hadrons, whereas the electron sample contains small 
contributions from Dalitz pairs and photon conversions, 
as well as some residual background from non-qg pro- 
cesses. We observe a clear difference in the normaliza- 
tion, not only in the relatively small event sample passing 
the neural-network selection, but also for the much larger 
sample available before the neural-network suppression. 
To correct for this difference, we apply additional scale 
factors to the simulated qq samples; they are different for 
electrons and muons. 

In addition to correcting the normalization, we also 
examine the shapes of the toeSi ^E, and q^ distribu- 
tions that are used to extract the signal yield. Since 
the size of the off-resonance data set is small, we study 
samples with a looser selection, namely we bypass the 
qq neural-network discrimination. The comparison of 
these gg-enriched samples reveals small differences be- 
tween data and simulation. We derive linear corrections 
from the bin-by-bin ratios and apply these corrections to 
the TTiEs, A£', and q^ distributions. Figures [T^ and [T51 



show a comparison of the rescaled and corrected qq MC 
samples with the off- resonance data for the AE", toes, 
and q^ distributions. Within the relatively large statisti- 
cal errors of the off-resonance data the simulation agrees 
well with the data. The uncertainties in the shape of the 
simulated distributions will be assessed as a systematic 
uncertainty. 



B. B ^f Xcf-u Enhanced Sample 

The overall dominant background source in this anal- 
ysis is B — )■ Xciv decays. Therefore it is important to 
verify that these decays are correctly simulated. This has 
been done in two ways, a) by relaxing the B — )■ Xciv 
suppression to obtain a charm-enhanced sample, and 
b) by reconstructing a specific decay mode, such as 
S" — > D*^£^iy, in the same way we reconstruct the signal 
decays, and comparing the kinematic distributions with 
MC simulations (see Section FVCp . 

We select a charm-enhanced sample by inverting the 
cut on the B — >■ Xcih' neural-network discriminator. Fig- 
ures [TJ] and [15] show the AE and tties distributions in the 
signal region and the side bands, as well as the g^ distri- 
bution in the signal region. All distributions show good 
agreement in shape; the absolute yields differ at a level 
that is expected; considering that the MC distributions 
have not been adjusted. 



C. B° -^ D'-£+u Control Sample 

To study the Monte Carlo simulation of the neutrino 
reconstruction employed in this analysis, we use a con- 
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FIG. 12: (color online) Comparison of off-resonance back- 
ground to B" — !> Tv^i'^v samples for data and MC distribu- 
tions. Top row: tties, center row: AE, and bottom row: q^ , 
separately for the electron (left column) and muon (right col- 
umn) samples. The shaded histograms indicate the true lep- 
tons, the hatched histograms indicate the fake leptons. The 
distributions are obtained from the full event selection, except 
for the qq neural-network discrimination. Linear corrections 
have been applied to the simulation. 



FIG. 13: (color online) Comparison of off-resonance back- 
ground to B" — > p^l'^u samples for data and MC distribu- 
tions. Top row: tties, center row: AE, and bottom row: q^ , 
separately for the electron (left column) and muon (right col- 
umn) samples. The shaded histograms indicate the true lep- 
tons, the hatched histograms indicate the fake leptons. The 
distributions are obtained from the full event selection, except 
for the qq neural-network discrimination. Linear corrections 
have been applied to the simulation. 



trol sample of exclusively reconstructed S" — > D*^ t^v 
decays with D*^ -^ D"tt~ and ZJ° -> K+n^ . Since 
the i?° -^ D*^i^v decay rate exceeds the rate for 
_B+ -^ p^l^v by a factor of about 30 (including the 
D*^ branching fraction), this control sample represents 
a high-statistics and high-purity sample of exclusive 
semileptonic decays. Except for the low-momentum pion 
(ttJ), this final state has the same number of tracks, and 
very similar kinematics, as the B'^ — >■ p'^l'^v signal decay. 
Furthermore, since about 50% of the B — )• Xctv back- 
ground in the B — )■ ■k(.v and B — >■ piv samples comes from 
5° -^ L»*"£+^ decays, this B° -> D*~£+i/ sample can 
provide important tests of the shapes of the distributions 
that are used to discriminate the B — >■ Xciv background 
from signal. 

Moreover, the distributions of the primary background 
suppression variables, in particular R2, L2, cos^by, 
^miss, and m^;gg/(2£^miss), are relatively insensitive to the 
specific semileptonic decay mode. Likewise, the resolu- 
tions for the fit variables AE, toeS; and q^ are dominated 
by the resolution of the reconstructed neutrino, and thus 
depend little on the decay mode under study. 

The reconstruction of the D*~ from its decay prod- 



ucts is straightforward. Except for the selection of the 
D*~ , we apply the same preselection as for the sig- 
nal charmless decays. We require the K~^n~ invariant 
mass to be within 17 MeV of the nominal D^ mass, 
and restrict the mass difference, Atoo* = mu^ — mu 
to 0.1432 < Amo' < 0.1478 GeV. The number of events 
in this data control sample exceeds the MC prediction by 
3.8 ± 1.7%, a result consistent with the uncertainties in 
the efficiency for the very-low-momentum charged pion 
from the D*^ -^ D^-k^ decay. We correct the MC yield 
and sequentially place requirements on the same seven 
variables we use in the neural networks to both the data 
and MC samples. We compare the step-by-step reduc- 
tion in the number of events; the largest difference is 
0.9 ± 0.7%, for the cut on cosOwt.- For all other critical 
requirements the agreement is better than 0.5% and one 
standard deviation. The remaining background is at the 
level of 10%. 

We have compared the MC-generated distributions for 
the control sample with the selected i?° -^ D*~l'^v data 
sample and find very good agreement for the basic event 
variables, i.e., the multiplicity of charged particles and 
photons, and the total charge per event, indicating that 



21 



(A 




B 


2.5 


a 




■a 


? 


■a 




c 

(0 


1.5 


O 


1 




0.5 







(3 




S 


1.2 




1.1 


CO 


1 




<^ 


0.9 
0.8 



xioi 



. xlO- 



5.1 5.15 



5.2 5.25 

m,s(GeV) 



■'5.1 



5.15 



5.2 5.25 
iHes (GeV) 




0.5 1 

AE (GeV) 




O 1-3 
E1.2 
^ 1.1 

5 1 

no.9 

°0.8 

0.7, 



5.2 5.25 

m,s(GeV) 



0.5 1 

AE (GeV) 



5.2 5.25 
rriEs (GeV) 



oi-3 

S1.2 

^ 1.1 
B 1 
rao.9 
°0.8 
0.7 



-0.5 



0.5 1 

AE (GeV) 



60 



40 



20 



Ol-3r 
S1.2 
^1.1 
S 1:* 
rao.9: 
°0.8 
0.7t— 



xlO' 


^ 


: 




-^ 



10 15 20 25 



10 15 20 25, 




xioi 



5.15 



5.2 5.25 
iHes (GeV) 



(00.25 



■°0.15 

(0 

O 0.1 
0.05 



Ol-3 
S1.2 
^ 1.1 

B 1 

no.9 

°0.8 

0.7, 




Ol-3 
S1.2 
^ 1.1 

B 1 

rao.9 

°0.8 

0.7 



-0.5 



x10° 






- 


• 













10 15 



20 



q2 (GeV^) 



0.5 1 

AE (GeV) 



T 


h,-'^_._. 


_fi_i_'__, 


1 


-0.5 


0.5 1 

AE (GeV) 



0.5 1 

AE (GeV) 



0.5 1 

AE (GeV) 



q^ (GeV^) 



15 



20 



FIG. 14: (color online) Comparison of data with MC distribu- 
tions for the charm-enhanced selection for the B° — >■ tt~£'^v 
sample. Top row: tties and AE for the signal bands, center 
row: rriES and AE for the side bands, and bottom row: Qcorr 
for the whole fit region. The bin-by-bin ratio of data over the 
sum of all MC contributions is given in the plots below each 
histogram. 



FIG. 15: (color online) Comparison of data with MC distribu- 
tions for the charm-enhanced selection for the B'^ — >■ p~ t^v 
sample. Top row: tties and AE for the signal bands, center 
row: rriES and AE for the side bands, and bottom row: Qcorr 
for the whole fit region. The bin-by-bin ratios of data over 
the sum of all MC contributions is given in the plots below 
each histogram. 



the efficiency losses are well reproduced by the simula- 
tion. The distributions of the topological event variables 
R2 and L2 match well. Figure [TBI shows the distributions 
of the variables critical for the neutrino reconstruction, 
m^igg/(2Smiss), cos6'by, and 6',„iss; they are also 



Pi 



well reproduced. 

Figure [T7] shows distributions of AE and ttjes for 
events in the signal region and in the side bands. Again, 
the agreement between data and the MC simulation is 



reasonable. 

We have also compared the q^ distributions of the sim- 
ulation and the data control sample and find good agree- 
ment for both the raw and the corrected spectra, as illus- 
trated in Figure [T51 After corrections, no events appear 
above the kinematic limit of 10.7 GeV . The q^^^^ resolu- 
tion function can be described by the sum of two Gaus- 
sian resolution functions, with widths of 0.27 GeV^ and 
0.67 GeV^, close to the values obtained for events in the 
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FIG. 16; (color online) Comparison of data and MC-simulated 
distributions for the B" — >■ D*^ t^v,D'^ — >■ A'+tt^ sample, af- 
ter selection criteria have been applied on all variables except 
the one presented, a) cos^by, b) Pmiss, c) 6miss, d) cos Owi, 
e) mniiss/2_Emiss, and f) the total charge per event Qtot. The 
background to the sample is indicated as a shaded (yellow) 
histogram. The combinatorial signal contribution is indicated 
as dashed histogram. 



fit region for the signal B 
respectively. 



■k(.u and B — > plv decays, 
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FIG. 17; (color online) Comparison of data and MC-simulated 
distributions for the 5° — >• D*~ (.^v,D° — >• K^-k~ sample, 
after all selection cuts have been applied, a) AiJ for events 
in the tties signal band, b) tties for events in the /S.E signal 
band, c) Ai? for events in the tties side band, and d) tties for 
events in the Ai? side bands. The background to the sample is 
indicated as a shaded (yellow) histogram. The combinatorial 
signal contribution is indicated as dashed histogram. 
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FIG. 18; (color online) Comparison of the data and MC simu- 
lation of q^ distributions for the B° ->• D*~£+u, D° ->• K+n~ 
sample after all selection criteria have been applied, a) the raw 
q^, and b) the corrected q^. The background to the sample is 
indicated as a shaded (yellow) histogram. The combinatorial 
signal contribution is indicated as dashed histogram. 



We determine the yields for the signal decay modes, 
B° -^ TT-£+i^, B+ -^ n°e+iy, B° -^ p~(-^v, and B+ -^ 
p^t^v, by performing a maximum-likelihood fit to the 
three-dimensional IS.E — tties — Q^ distributions for the 
four selected data samples corresponding to the four ex- 
clusive decay modes. The fit technique employed in this 
analysis is an extended binned maximum-likelihood fit 
that accounts for the statistical fiuctuations not only of 
the data samples but also of the MC samples by allow- 
ing the MC-simulated distributions to fluctuate in each 
bin according to the statistical uncertainty given by the 
number of events in the bin. This method was introduced 
by Barlow and Beeston [4g |. 

The parameters of the fit are scale factors for the signal 
and background yields of the four selected event samples. 
We use the following nomenclature for the fit parame- 



ters: p™'"'''^''^ where the superscript source denotes the 
fit source (signal or background type) and the subscript 
j labels the q^^^^ bin (if no subscript is given, the same fit 
parameter is used across all q^ bins). Predictions for the 
shape of the Ai?— mES distributions are taken from simu- 
lation of both signal and the various background sources, 
separately for each bin in q^. The branching fractions 
for the four signal decays are obtained by multiplying 
the fitted values of the scale factors with the branching 
fractions that are implemented in the MC simulation. 

The choice of a two-dimensional distribution in /S.E 
and TTiES is mandated because the two variables are corre- 
lated for both signal, in particular the combinatorial sig- 
nal events, and for some of the background sources. Since 
it would be difficult to determine reliable analytic expres- 
sions for these two-dimensional distributions, a binned 
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maximum-likelihood method is used, with the bin sizes 
chosen to obtain a good signal and background separa- 
tion while retaining adequate statistics in all bins. The 
bin sizes are small in the region where most of the sig- 
nal is located and larger in the side bands. There are 47 
AE — niES bins for each bin in (/con- Figure [T9l shows the 
AE — Toes distribution for signal events and the binning 
used in the fit. As mentioned in Section flVB 21 for the 
two B — )• TT£iy samples the g^ range < q^^^.^ < 26.4 GcV^ 
is divided into six bins, and for the two B — > pliy sam- 
ples the range < q'^^^^ < 20.3 GeV^ is divided into three 
bins. 
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FIG. 19: Distribution of AE vs twes for true B° -i> tt^^+i/ 
signal events (left) and definition of bins in the AE — tues 
plane used in the fit for all samples (right). 



B. Fit Method 

Since the MC samples available to create the probabil- 
ity density functions (PDFs) for the individual sources 
that are input to the fits are rather limited in size, it is 
necessary to take into account the statistical uncertain- 
ties, given by the number of events generated for each bin. 
For this reason we have adopted a generalized binned 
maximum-likelihood fit method. The MC samples that 
are used to define the PDFs are to a good approximation 
statistically independent of those used to train the neural 
networks for background suppression, since for the latter 
relatively small subsamples of the full MC samples have 
been used. 

As mentioned above, the data are divided into n bins 
in a three-dimensional array in AE — tties — 'Zcorr- 

If di is the number of selected events in bin i for a 
given single data sample corresponding to candidates for 
a specific decay mode, and Uji is the number of MC events 
from source j in this bin, then 



Nr 



-Ed^, 



N, 






ijz, 



(24) 



where Nd is the total number of events in the data sam- 
ple, and Nj is the total number in the MC sample for 
source j. We assume that there are m different MC- 
generated source distributions that add up to describe 
the data. The predicted number of events in each bin 



fi{Pj) can be written in terms of the strength of the in- 
dividual contributions Pj (j = 1, ..,to) as 



ft = NdY^^ PjWj^aj^/Nj = VpjWjzaj,, 



/ 



/ , 1 

3=1 



(25) 



with Pj ~ NdPj/Nj. In each bin, the weights Wji ac- 
count for the relative normalization of the samples and 
various other corrections. 

Since the MC samples are limited in size, the generated 
numbers of events aji have statistical fluctuations relative 
to the value Aji expected for infinite statistics, and thus 
the more correct prediction for each bin is 






PjWj.A^ 



ji- 



(26) 



If we assume Poisson statistics for both the data and MC 
samples, the total likelihood function C is the combined 
probability for the observed di and aji [49| . 

n n m 

In £ = Y,id^ In f^ - h) + E E("j^ ^^^ ^1^ - ^J')' (^7) 

i=l i=l j=l 

The first sum has the usual form associated with the un- 
certainty of the data and the second term refers to the 
MC statistics and is not dependent on data. There are 
(n -|- 1) X TO, unknown parameters that need to be de- 
termined: the m relative normalization factors pj , which 
are of interest to the signal extraction, and n x to values 
Aji. 

The problem can be significantly simplified. The nxm 
quantities Aji can be determined by solving n simulta- 
neous equations for Aji of the form 



2_^PjWjiAji 



PjWj 



^^ 1 



Pj^ji'^i 



(28) 



with Aji ~ aji/ [1 



PjWjiU) and U = l-di/ Ji (for di = 
we define t; = 1). At every step in the minimization of 
— 21n£ these n independent equations need to be solved. 
This procedure results not only in the determination of 
the parameters pj , but also in improved estimates for the 
various contributions Aji in each bin. 

For fits to the individual data samples corresponding to 
the four signal decay modes, there is a specific likelihood 
function (Eq. [27)) . To perform a simultaneous fit to all 
four data samples the log-likelihood function is the sum 
of the individual ones. Some of the parameters pj may 
be shared among the four likelihood functions. 



ln/: = ^ln£„ = ^5I(^'"ln/''^-^*) (29) 

4 71 m 
h=l i=l j=l 
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C. Fit Parameters and Inputs 

The fits can be performed separately for each of the 
four data samples or combined for all four data sam- 
ples, and where possible, with common fit parameters. 
The nominal fit in this analysis is a simultaneous fit of 
all four data samples: 5° -^ n^l+i/, B+ -^ 7r°^+i^, 



B^ 



-£+v. and B+ 



p^l^v. A signal decay in 



one data sample may contribute to the background in 
another sample, and therefore these sources share a com- 
mon fit parameter. For example, the scale factor for the 
B^ -^ p^l'^v signal in the B^ — > p^l^v sample is also 
applied in the B^ — >■ 7r~^+zy sample, where it represents 
cross-feed background. 

We impose isospin invariance for the signal decay 
modes, 



r(B° ^ p-i+v) = 2T{B^ 



p"e+i^) . (30) 



The yields of the true and combinatorial signal decays as 
well as isospin-conjugate decays are related to the same 
branching fraction and therefore share the same fit pa- 
rameter. The B — >■ X^lv background, which contains 
exclusive and non-resonant decays, is scaled by two pa- 
rameters, one for low and intermediate q^ {q^ < 20 GeV^) 
and one for high q"^ {q^ > 20GeV^), for the fits to the 
B — )• nlv samples. Because of the large correlation 
between B — !• X^tv background and B — )■ p(,v signal 
(> 90% for both B -^ ptv modes), we rely on MC sim- 
ulation for the B -^ Xutv background and keep it fixed 
in the fits to the B — >• p(.v samples. The BB background 
is split into two sources. Among the B — ?► XJm de- 
cays, we treat the dominant decay mode, B — > D*iiy, 
as a separate source and combine the other semileptonic 
decays {B — >■ Dli/, B -> D^*\mT)(.v) and the remain- 
ing (or "other") BB background (secondary leptons and 
fake leptons) into a single source. The continuum qq 
background sources containing true and fake leptons arc 
combined into one fit source and scaled by a single fit 
parameter. 

The complete list of fit sources and corresponding fit 
parameters is given in Tabic FvTl The tt o p cross feed is a 
free fit parameter in the four-mode fit; for one-mode fits, 
it is fixed to the value obtained from the four- mode fit. 
In the four- mode fit, all background sources that arc not 
fixed arc fit separately for each signal mode, since the 
different hadrons of the signal decays lead to different 
combinatorial backgrounds. 



D. Fit Results 

The fits are performed both separately and simulta- 
neously for the four signal decay modes, B^ — > 7r~^+i^, 
B+ -^ TT°i+iy, B° ^ p-i+iy, and B+ -^ p"e+i^. Fig- 
ures [5DH23] show projections of the fitted AE — tties dis- 
tributions in the signal bands for these decays, separately 



TABLE VI: List of fit parameters representing scale factors 
for the difi'erent signal samples and background sources. Pa- 
rameters with index j are free parameters in the fit, one for 
each q^ bin j. The n <r^ p crossfeed parameter is free only 
in the four- mode fit; for one-mode fits, it is fixed to the val- 
ues obtained from the four-mode fit. There are independent 
scale factors for qq background, B — >■ D*£u decays and for all 
other background sources from BB events for all four signal 
modes (subscripts 7r^,7r°,p^,p''). For the B — >■ tt£u decays, 
the B — >■ Xu£u background is fit in two q^ intervals (index 



fc = 1, 2); for the B — >• ptv decays it is fixed. 




Source / sample B — )> tv£u 


B -^ plu 


Q;n-,ie,1 «'^'^'' 


ptv 



pf 
pf 

plv 

Pj 

Pn±,k' PnO,k 

B ^ D*(.v background p^l'^" , p^o'*" 



Combinatorial signal 
Isospin-conjugate signal 
Cross feed tt -f-)- p 
B — >■ Xulv background 



Pj 

ptv 

Pj 

plv 

Pj 
pf 

fixed 

„D*lv „D*li/ 
Po± . PoO 



Other BB background pfi^"""^^, p°'o'"="-^^ Ppi""^^, P^. 



therBB 





qq background pl^ , p™ p'i , P^o 



for each bin in Qcorr- ^^ a measure of the goodness-of-fit 
we use x^ P^r degree of freedom; all fits have values in 
the range 1.05-1.11 (for details see Table EIIl)- 

The scale factors for the signal contributions, which 
are determined by the fits, can be translated to numbers 
of background-subtracted signal events for the four signal 
decays. These signal yields are listed in Table FVlIl with 
errors that are a combination of the statistical uncertain- 
ties of the data and MC samples and the uncertainties 
of the fitted yields of the various backgrounds. For each 
signal decay mode, the table specifies the number of true 
and combinatorial signal decays. Their relative fraction 
is taken from simulation. This fraction is larger for de- 
cays with a 7r° in the final state. For all signal modes, 
the fraction of combinatorial signal events is small at low 
q^, increases with q^, and at the highest q^ it is similar 
to or exceeds the one of true signal decays. This leads to 
larger errors in the measurement of q'^, toes and AE. 

In Tablc lXVlIl in Appendix IXI BI thc correlation matrix 
of the four-mode fit is presented. We observe correlations 
of about 40 — 60% between the qq and the other BB 
backgrounds and between the B — >■ D*ii' and the other 
BB backgrounds for all signal modes. For B — )- nii^, 
the correlation between the B — >■ Xuii' background and 
the signal at high q"^ is also sizable (~ 60%). For B ^■ 
plv, this correlation is larger than 90%, which is why we 
choose to fix the B — > Xuiiy background normalization 
for these two samples. As a test, we let the B —>■ Xuiv 
background normalization in the B — ^ piv modes vary as 
free parameter in the four- mode fit. This results in a i? — >■ 
Xuiv contribution that is lower by a factor of 0.85 ±0.15 
for S" -> p-l+u and 0.90 ± 0.14 for B+ -^ p^l+u and 
an increase of the B — !• piv signal yields by 10% in the 
first two q^ bins and by 15% in the last q^ bin. These 
changes are covered by the systematic uncertainties due 
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to the B — >■ Xuiiy background stated in Section IVlII 

To cross-check the results of the nominal four-mode 
fit, we also perform fits for each signal mode separately. 
The contributions from the other signal decay modes are 
fixed to the result obtained from the four-mode fit. Since 
the shape of the tt o p cross-feed contribution is very 
similar to the other B — > Xutv background, we fix its 
normalization to the one obtained from the four-mode 
fit. A comparison of the results of the one-mode fits with 
the combined four-mode fit shows agreement within the 
fit errors of the B^ — >■ 7r~^+z/ and B'^ — >■ 7r°£+i/ modes 
and the i?" -^ p^i^iy and _B+ -^ p'^t^v modes in all q^ 
bins. 

The partial branching fractions for the different (j'^orr 
bins arc derived as the products of the fitted signal scale 
factors and the signal decay branching fractions used in 
the simulation. The total branching fraction integrated 
over the entire q^ range and its error are calculated as 
the sum of all partial branching fractions, taking into 
account the correlations of the fitted yields in different 
q^ bins. The branching fraction for B^ decays, B^-^ ^, is 
related to the fitted signal yields, -/V°j j, in the following 
way. 



signal 



signal 



4 X e 



^0 f, 

signal •' ^ 



OO^BB 



(31) 



where /oo == 0.484 ± 0.006 [4l| denotes the fraction of 
B^B'^ events produced in T(45) decays and Cgignai is the 
total signal efficiency (averaged over the electron and 
muon samples) as predicted by the MC simulation. The 
factor of four accounts for the fact that each event con- 
tains two B mesons, and that the branching fraction is 
quoted for a single charged lepton, not for the sum of the 
decays to electrons or muons. The branching fraction 
results are presented in Section IVIIII 



E. Fit Validation and Consistency 

The fit procedure is validated several ways. First of all, 
the implementation of the Barlow-Becston fit technique 
allowing statistical fluctuations of the MC distributions 
to be incorporated is checked by verifying the consistency 
of the fit variations with the statistical error of the input 
distributions. Secondly, a large number of simulated ex- 
periments are generated based on random samples drawn 
from the three-dimensional histograms used in the stan- 
dard fit. Specifically, we create 500 sets of distributions 
by fluctuating each simulated source distribution bin-by- 
bin using Poisson statistics. For each of the sets, we 
add the source distributions to make up to total distri- 
bution that corresponds to the data distribution ("toy 
data"), which are then fitted by the standard procedure. 
In addition, we create independent fiuctuations for the 
distributions that make up the source PDFs for the fit, 
in the same way as for the toy data described above. For 
a compilation of these 500 "toy experiments" , we study 



the distributions of the deviation of the fit result from the 
input value divided by the fit error. These distributions 
show no significant bias for any of the free parameters 
and confirm that the errors are correctly estimated. 

Additional fits are performed to check the consistency 
of the data. For instance, the data samples are divided 
into subsamples, i.e., the electron sample separated from 
the muon sample or the data separated into different run 
periods. These subsamples are fitted separately; the re- 
sults agree within the statistical uncertainties. 



VII. SYSTEMATIC UNCERTAINTIES 

Many sources of systematic uncertainties have been as- 
sessed for the measurement of the exclusive branching 
fractions as a function of q^. Since this analysis does not 
depend only on the reconstruction of the charged lepton 
and hadron from the signal decay mode, but also on the 
measurement of all remaining tracks and photons in the 
event, the uncertainties in the detection efficiencies of all 
particles as well as the uncertainties in the background 
yields and shapes enter into the systematic errors. 

Tables IVIIII and IIXI summarize the systematic uncer- 
tainties for B — > Triiy and B — ?► piv for the four-mode 
fit. In Appendix IXI Al the systematic error tables for the 
one-mode fits are presented. The individual sources are, 
to a good approximation, uncorrelated and can therefore 
be added in quadrature to obtain the total systematic 
errors for each decay mode. In the following, we discuss 
the assessment of the systematic uncertainties in detail. 

For the estimation of the systematic errors of the fitted 
branching fractions, we compare the differential branch- 
ing fractions obtained from the nominal fit with results 
obtained after changes to the MC simulation that refiect 
the uncertainty in the parameters that impact the detec- 
tor efficiency and resolution or the simulation of signal 
and background processes. For instance, wc vary the 
tracking efficiency, reprocess the MC samples, reapply 
the fit to the data, and take the difference compared to 
the results obtained with the nominal MC simulation as 
an estimate of the systematic error. The sources of sys- 
tematic errors are not identical for all four signal decay 
modes, and the size of their impact on the event yields 
depends on the sample composition and q^. 



A. Detector Effects 

Uncertainties in the reconstruction efficiencies for 
charged and neutral particles and in the rate of tracks 
and photons from beam background, fake tracks, failures 
in the matching of EMC clusters to charged tracks, and 
showers split off from hadronic interactions, undetected 
Kl, and additional neutrinos, all contribute to the qual- 
ity of the neutrino reconstruction and impact the vari- 
ables that are used in the preselection and the neural 
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FIG. 20: (color online) rriES and AE distributions in each q^ bin for B" — >■ tt I'^iy after the fit. The distributions are shown in 
the AE and tties signal bands, respectively. Legend: see Figure [S] 



networks. For all these effects the uncertainties in the ef- 
ficiencies and resolution have been derived independently 
from comparisons of data and MC simulation for selected 
control samples. 



1. Track, Photon, and Neutral-Pion Reconstruction 

We evaluate the impact of uncertainties in the tracking 
efficiency by randomly eliminating tracks with a proba- 



bility that is given by the uncertainty ranging from 0.25% 
to 0.5% per track, as measured with data control sam- 
ples. 

Similarly, we evaluate the uncertainty due to pho- 
ton efficiency by eliminating photons at random with 
an energy-dependent probability, ranging from 0.7% per 
photon above 1 GeV to 1.8% at lower energies. This es- 
timate includes the uncertainty in the 7r° efficiency for 
signal decays with a tt", since photons originating from 
the signal hadron are also eliminated. 
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FIG. 21: (color online) tties and AE distributions in each q^ bin for B^ — >■ Tr"^^!/ after the fit. The distributions are shown in 
the AE and tties signal bands, respectively. Legend: see Figure [S] 



2. Lepton Identification 



3. A'? Production and Interactions 



The average uncertainties in the identification of elec- 
trons and muons have been assessed to be 1.4% and 3%, 
respectively. The uncertainty in the misidcntification of 
hadrons as electrons or muons is about 15%. 



Events containing a A'^ have a significant impact on 
the neutrino reconstruction, because only a small frac- 
tion of the K^ energy is deposited in the electromagnetic 
calorimeter. Based on detailed studies of data control 
samples of I?° — > K'^tt^tt~ decays and inclusive K'g sam- 
ples in data and MC, corrections to the efficiency, shower 
deposition and the production rates have been derived 
and applied to the simulation as a function of the K'^ mo- 
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FIG. 22: (color online) tties and AE distributions in each q^ bin for _B" 
the AE and tties signal bands, respectively. Legend: see Figure O 



v after the fit. The distributions are shown in 
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mentum and angles (see Section |III|) . To determine the 
systematic uncertainties in the MC simulations we vary 
the scale factors within their statistical and systematic 
uncertainties. The average uncertainty of the energy de- 
position in the EMC due to K^ interactions is estimated 
to be 7.5%. Above 0.7 GeV, the K^ detection efficiency 



is well reproduced by the simulation, with an estimated 
average uncertainty of 2%. At lower momenta, the simu- 
lation is corrected to match the data, and the uncertainty 
increases to 25% below 0.4 GeV. 

The production rates for K^ in data and MC agree 
within errors, except for momenta below 0.4 GcV where 
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TABLE VII: Results of fits performed separately for each of the four signal decay modes and simultaneously for all four decay 
modes in bins of q^ : x^ per degree of freedom, signal yields for true signal decays, 7V°'^ , and combinatorial signal, TV'^"™'^ . The 
stated errors are the fit errors. 



q^ range (GeV^) 



X^ndf 0-4 



4-8 



8-12 



12-16 



16-20 



> 20 < g^ < 26.4 






N„± = N^'^ + N'^T*' 



259/268 701 ± 64 1950 ± 104 
IzbO.l 4 ±0.2 
702 ± 64 1954 ± 104 



1552 ±113 1184 ±104 

9 ± 1 30 ± 3 
1561 ±113 1214 ±104 



732 ± 80 541 ± 68 

77 ±8 401 ±51 

809 ± 80 942 ± 85 



6660 ± 2 78 

521 ± 22 
7181 ± 279 






237/268 315 ±42 576 ± 54 
5 ±0.7 13 ±1 
320 ± 42 589 ± 54 



904 ± 107 471 ± 68 

36 ± 4 36 ± 5 

940 ± 107 507 ± 68 



414 ± 83 159 ± 55 
119 ±24 397 ±137 
533 ± 86 556 ± 148 



2840 ± 203 

606 ± 43 
3446 ± 208 



N^± + N^o 4-mode 


799/819 1012 ± 76 2535 ± 128 2485 ± 157 1729 ± 130 1291 ± 125 1552 ± 180 10604 ± 376 


g= range (GeV^) 


0-8 8-16 >16 < q^ < 20.3 






Np± = n:i' + iv„± 



147/131 



comb 

pi 



237 ± 56 
56 ±13 
293 ± 57 



459 ± 44 
287 ± 27 
746 ± 52 



170 ±17 
368 ± 38 
538 ± 42 



866 ± 101 
711 ± 82 
1577 ± 130 



\TComb 



162/131 



253 ± 63 

31 ±8 

284 ± 64 



856 ± 74 
267 ± 23 
1123 ± 77 



294 ± 29 
270 ± 27 
564 ± 40 



1403 ± 143 

567 ± 58 

1970 ± 154 



Np± + NpO 4-mode 



799/819 



471 ± 101 



1754 ± 120 



1105 ±: 



3332 ± 286 



the data spectrum is low by 22 ± 7% compared to the 
MC simulation and a correction is applied. To assess the 
impact of the uncertainty of the correction procedure, 
the size of the correction is varied by its estimated un- 
certainty. 



the difference between LCSR and ISGW2 as systematic 
uncertainty, but this difference is comparable to the es- 
timate we obtain from the uncertainties in the LCSR 
calculation. 



B. Simulation of Signal and Background 

1. Signal Form Factors 

To assess the impact of the form-factor (FF) uncer- 
tainty on the shape of the simulated signal distributions, 
we vary the B — ^ tt form factor within the uncertainty 
of the previous BaBAR measurement Q and the B ^ p 
form factors within the uncertainties of the LCSR calcu- 
lation assessed by Ball and Zwicky [131 ■ ^^^ the latter we 
assume uncertainties on the form factors Ai , A2 and V of 
10% at g2 = 0. They rise hnearly to 13% at q'^ = 14 GeV^ 
and are extrapolated up to the kinematic endpoint. We 
add the uncertainties due to the three form factors in 
quadrature. For B — >■ tt£v, the form-factor uncertainty 
is small, since we extract the signal in six bins of q^. 
In contrast, for B — > pliy the form- factor uncertainty is 
one of the dominant sources of systematic error. This 
is partly due to the stricter requirement on the lepton 
momentum, Pf > 1.8 GcV, which is imposed to suppress 
the large B — ^ Xclv background. We refrain from using 



8. _B — >■ Xutv Background 

The B — > Xuiv background contribution is composed 
of the sum of exclusive decays, B^ — > ojI'^v, B^ — > rjl^i', 
and B~^ — )■ rj'i'^v decays, and the remaining resonant and 
non-resonant B — )■ Xuii' decays that make up the total 
B — )• Xuii' branching fraction. We estimate the total er- 
ror of the B — > Xuiv background composition by repeat- 
ing the fit with branching fractions for various exclusive 
and non-resonant decays varied independently within 
their current measurement errors. The uncertainty of the 
branching fraction for non- resonant decays is dominant; 
it is equal to the error on the total B — > X^f^i/ branching 
fraction, B{B -^ X^tv) = (2.33 ± 0.22) x lO'^ [4lj. 

In addition, the analysis is sensitive to the mass and 
composition of the charmless hadronic states. We assess 
the uncertainty of the predictions by varying the QCD 
parameters that define the mass, the lepton spectrum, 
and the q^ distributions predicted by calculations [3Q| 
based on HQE. We vary the shape-function (SF) param- 
eters mil and ui within the uncertainties (error ellipse) 
given in Ref. [38|. 



TABLE VIII: Systematic errors in % for Z3(B" -^ n-l+v) 
from the four-mode fit for bins in q and the total q range. 
The total errors are derived from the individual contributions 
taking into account the complete covariance matrix. 
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TABLE IX: Systematic errors in % for B{B° -^ p'l+u) from 
the four-mode fit for three bins in q and the total q range. 
The total errors are derived from the individual contributions 
taking into account the complete covariance matrix. 





B 


' — >■ 


TTliy 










B -) 


. plv 








q-' range (GeV^) 


0-4 4-8 8-12 


12-16 


16-20 


>20 0-26.4 


q'' range (GeV") 
Track efficiency 


0-8 
3.2 


8-16 
2.9 


>16 
0.3 


0-20.3 


Track efficiency 


3.4 


1.5 


2.3 


0.1 


1.5 


2.8 


1.9 


2.5 


Photon efficiency 


0.1 


1.4 


1.0 


4.6 


2.8 


0.3 


1.8 


Photon efficiency 
Lepton Identification 
Kl efficiency 


2.6 
5.7 
10.3 


2.0 
3.0 
1.2 


2.6 
4.0 
4.9 


2.4 


Lepton identification 


3.8 


1.6 


1.9 


1.8 


1.9 


3.0 


1.8 


3.4 


Kl efficiency 


1.0 


0.1 


0.5 


4.5 


0.4 


2.0 


1.4 


4.8 


Kl shower energy 


0.1 


0.1 


0.1 


0.8 


0.9 


3.8 


0.7 


Kl shower energy 


1.6 


0.8 


1.0 


1.1 


Kl spectrum 


1.6 


1.9 


2.2 


3.1 


4.4 


2.3 


2.5 


Kl spectrum 
B -^ Tvlu FF /+ 


4.2 
0.1 


6.1 
0.1 


7.0 
0.7 


5.7 


B -^ -kIuFF f+ 


0.5 


0.5 


0.5 


0.6 


1.0 


1.0 


0.6 


0.2 


B -> pli^FFAi 


1.7 


1.2 


3.4 


2.0 


0.1 


1.6 


1.7 


B ->• plf FF Ai 


10.7 


6.6 


4.5 


7.5 


B ->• p£uFFA2 


1.3 


0.8 


2.6 


1.0 


0.1 


0.4 


1.1 


B -^ plv FF A2 


8.5 


3.8 


0.8 


4.7 


B -^ pluFFV 


0.2 


0.3 


0.9 


0.7 


0.1 


0.5 


0.5 


B -> plu FF V 
B{B+ -^ Lol+u) 


3.4 
0.7 


3.0 
0.7 


3.6 
3.4 


3.2 


B{B+ -^ uj£.+u) 


0.1 


0.1 


0.1 


0.2 


0.3 


1.5 


0.2 


1.2 


B{B+ -^ rji+u) 


0.1 


0.1 


0.2 


0.2 


0.2 


0.5 


0.2 


B{B+ -> r]e+u) 


0.8 


0.1 


0.6 


0.4 


B{B+ -^ rjtv) 


0.1 


0.1 


0.1 


0.1 


0.1 


0.3 


0.1 


B{B+ -^ r]'£+u) 


0.8 


0.5 


1.2 


0.7 


B{B -^ Xulu) 


0.2 


0.1 


0.1 


0.1 


1.1 


1.6 


0.4 


B{B -^ Xulv) 


7.4 


7.3 


10.6 


8.0 


B — >• X^iv SF param. 


0.4 


0.1 


0.2 


0.2 


0.5 


4.2 


0.7 


B — >• X^Hv SF param. 
B -^ Dtv FF p^ 


11.9 
0.9 


7.6 
0.2 


12.8 
0.1 


10.0 


B -^ Dlv FF pi, 


0.2 


0.1 


0.5 


0.3 


0.2 


0.7 


0.3 


0.4 


B -> D*£iy FF i?i 


0.1 


0.4 


0.8 


0.6 


0.3 


0.6 


0.5 


B -> D*lu¥¥ Ri 


0.7 


0.1 


0.3 


0.3 


B -^ D*£iy FF R2 


0.5 


0.2 


0.1 


0.2 


0.1 


0.4 


0.2 


B -^ D*£uFF R2 


1.7 


0.1 


0.2 


0.6 


B -^ D*£v FF pI,. 


0.7 


0.2 


0.6 


0.8 


0.4 


1.1 


0.6 


B -> D*euFF p%. 


2.0 


0.2 


0.1 


0.7 


B{B -^ Dtv) 


0.2 


0.2 


0.3 


0.4 


0.5 


0.5 


0.3 


B{B ^ Dlu) 


1.6 


0.3 


0.1 


0.7 


B{B ^ D*iv) 


0.4 


0.1 


0.3 


0.3 


0.3 


0.7 


0.3 


B{B ^ D-lu) 


0.5 


0.1 


0.3 


0.3 


B{B -> D**^i/)narrow 


0.4 


0.1 


0.1 


0.3 


0.1 


0.5 


0.2 


B{B -> D"^//)narrow 


1.3 


0.1 


0.1 


0.5 


B{B -^ D"^//)teoad 


0.1 


0.1 


0.1 


0.5 


0.1 


0.2 


0.2 


B{B -^ D"^//)teoad 


0.7 


0.1 


0.1 


0.3 


Secondary leptons 


0.5 


0.2 


0.3 


0.2 


0.2 


0.7 


0.3 


Secondary leptons 

Continuum 

Bremsstrahlung 


1.5 
8.9 
0.9 


0.1 
3.8 
0.1 


0.1 
5.0 
0.2 


0.5 


Continuum 


5.3 


1.0 


2.6 


1.8 


3.1 


6.1 


2.0 


4.0 


Bremsstrahlung 


0.3 


0.1 


0.1 


0.1 


0.1 


0.4 


0.2 


0.4 


Radiative corrections 


0.5 


0.1 


0.1 


0.2 


0.2 


0.6 


0.3 


Radiative corrections 
Nbb 


1.3 
2.7 


0.1 
2.0 


0.7 
2.5 


0.6 


Nbb 


1.2 


1.0 


1.2 


1.2 


1.1 


1.6 


1.2 


2.3 


B lifetimes 


0.3 


0.3 


0.3 


0.3 


0.3 


0.7 


0.3 


B lifetimes 


1.5 


0.4 


0.4 


0.7 


/±//oo 


1.0 


0.4 


0.8 


0.8 


0.5 


1.3 


0.8 


/±//oo 

Total syst. error 


1.2 
26.1 


0.1 
16.1 


0.1 
21.3 


0.4 


Total syst. error 


8.2 


3.9 


6.7 


8.3 


6.9 


10.6 


5.0 


15.7 



For the two B — > pii/ samples, the B — )• Xuli^ back- 
ground is large compared to the signal and very difficult 
to separate. Consequently, the fit shows very high corre- 
lations between the fitted yields for signal and this back- 
ground. We therefore choose to fix the background yields 
and shapes to those provided by the simulation, and ac- 
count for the uncertainty by assessing the sensitivity of 
the fitted signal yield to variations of the B — )■ X^liy 
branching fraction and the shapes of the background dis- 
tributions, corresponding to the estimated error of the 
shape-function parameters. The resulting estimated er- 
rors are the two dominant contributions to the system- 
atic errors of the B — )• plv partial and total branching 
fractions. 



3. -B — > Xcliy Background 

The systematic error related to the shapes of the 
B — > Xcii" background distributions is dominated by the 



uncertainties in the branching fractions and form factors 
for the various semileptonic decays. We vary the compo- 
sition of the B — >■ Xciiy background based on a compila- 
tion of the individual branching fractions of -B — ?> D£i', 
B -^ 0*11/ and B -^ 0**11/ (narrow and broad D** 
states) decays within the ranges given by their errors, see 
Table m Since we scaled up the four B — > D**iv branch- 
ing fractions to take into account the unknown D** par- 
tial branching fractions, the errors were increased by a 
factor of three relative to the published values. 

To evaluate the effect of uncertainties in the form- 
factor parameters for the dominant B — >■ D*^i/ compo- 
nent, we repeat the fit with ztlcr variations in each of 
the three form-factor parameters, p^,,Ri and i?2- The 
impact of the form factor for the B -^ Dtv background 
is evaluated by varying the parameter p|) within its un- 
certainty. 
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4- Continuum Background 



In Section FVAl we have described the correction of the 
simulated shapes of the mES, Ai^, and q^ distributions 
for the continuum using hnear functions derived from 
comparison with off-resonance data. The uncertainties 
of the fitted slopes of these correction functions are used 
to evaluate the errors due to modeling of the shape of 
the continuum background distributions. They represent 
a sizable contribution to the systematic error, which is 
mainly due to the low statistics of the off-resonance data 
sample. 



C. Other Systematic Uncertainties 

1. Final-State Radiation and Bremsstrahlung 

The kinematics of the signal decays are corrected 
for radiative effects such as final-state radiation and 
bremsstrahlung in detector material. 

In the MC simulation, final-state radiation (FSR) is 
modeled using PHOTOS [H, which is based on 0{a) 
calculations but includes multiple-photon emission from 
the electron. We have studied the effects of FSR on the 
q^ dependence of the measured signal and background 
yields by comparing events generated with and with- 
out PHOTOS. The observed change is largest, up to 
5%, for electron momenta of about 0.6 GeV (i.e. well 
below our cut-off at 1 GeV for B — >■ nlv and 1.8 GcV 
for B — > p£v). Comparisons of the PHOTOS simula- 
tion with semi-analytical calculations [50| show excellent 
agreement. Allowing for the fact that non- leading terms 
from possible electromagnetic corrections to the strong 
interactions of the quarks in the initial and final state 
have not been calculated to any precision [51j . we adopt 
an uncertainty in the PHOTOS calculations of 20%. 

The uncertainty of the bremsstrahlung correction is 
determined by the uncertainty of the amount of detector 
material in the inner detector. Wc have adopted as the 
systematic uncertainty due to bremsstrahlung the impact 
of a change in the thickness of the detector material by 
0.14% radiation lengths, the estimated uncertainty in the 
thickness of inner detector and the beam vacuum pipe. 
As for final-state radiation, the uncertainty in the ef- 
fective radiator thickness impacts primarily the electron 
spectrum. 

The uncertainties due to final-state radiation and 
bremsstrahlung combined amount to far less than 1% for 
most of the g^ range. 



2. Number of BE Events 

The determination of the on-resonance luminosity and 
the number of BB events is described in detail else- 
where (52| . The uncertainty of the total number of BB 
pairs is estimated to be 1.1%. 



At the T[AS) resonance, the fraction of B^B'^ events 
is measured to be /oo = 0.484 ± 0.006, with the ratio 
/+-//00 = 1.065 ± 0.026 [m. This error impacts the 
branching-ratio measurements by 0.8%. 



3. B° and B+ Lifetimes 

Since we combine fits to decays of charged and neutral 
B mesons and make use of isospin relations, the i?-mcson 
lifetimes enter into the four-mode fit. We use the PDG 
Q value for the B lifetime, tq = 1.530±0.009 ps , and the 
lifetime ratio, tj^Jtq = 1.071 ±0.009. These uncertainties 
lead to a systematic error of 0.3% for B — > iriiy and 0.7% 
for B — > pii' decays. 



VIII. RESULTS 

Based on the signal yields obtained in the four-mode 
fit, integrated over the full g^ range (see Table IVIip . 
we derive the following total branching fractions, con- 
strained by the isospin relations stated in Eqs. I30[ 

B{B° ^TT-£+iy) = (1.41 ± 0.05 ± 0.07) X 10"" , 
BiB° ^ p-£+iy) = (1.75 ± 0.15 ± 0.27) X 10"" . 

Here and in the following, the first error reflects the sta- 
tistical (fit) error and the second the estimated system- 
atic error. The total branching fractions obtained from 
the single-mode fits for the charged and neutral B — >■ Triiy 
samples are 

B(S° -> 7r-£+j/) = (1.44 ± 0.06 ± 0.07) X 10"" , 
B{B+ ^n°e+iy) x2— = (1.40 ± 0.10 ± 0.11) x 10"" . 

T+ 

For the charged and neutral B -> pii/ samples, we obtain 

B{B° ^ p-e+iy) = (1.98 ± 0.21 ± 0.38) X 10"" , 
B{B+ ^ p^tv)-x2— = (1.87 ± 0.19 ± 0.32) X 10"" . 



The single- mode fits result in higher values for B{B'^ — >■ 
p^l'^v) and B{B^ — )■ p^i^v) than the average branching 
fraction obtained from the four-mode fit. This may be 
explained by different treatments of the isospin-conjugate 
signal and the tt <-> p cross feed in the single- and four- 
mode fits. In contrast to the four-mode fit, the isospin- 
conjugate signal contribution in the single-mode fits is 
not constrained by the isospin-conjugate mode. In ad- 
dition, the four-mode fit uses the same fit parameter 
for the signal and the cross feed from the signal mode 
into other modes, which leads to a slight decrease in the 
B — > piv branching fraction compared to the single-mode 
fits. Since the p ^- n cross feed is significantly larger than 
the TT — )• p cross feed, the effect on the B -^ plv results 
is larger than for B — )■ ttIv. 
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Both the B -^ nil' and the B — >■ piv results are con- 
sistent within errors with the isospin relations, 



B(B+ -^ TT^i+v) X 2^ 

B{B" -^ p-(+v) 
B{B+ -^ p"£+u) X 2^ 



1.03 ±0.09 ±0.06 



1.06 ±0.16 ±0.08 



By extracting the signal in several q^ bins we also mea- 
sure the q^ spectra of i? — >■ irliy and B -^ ptv decays. 
These spectra need to be corrected for effects such as 
detector resolution, bremsstrahlung, and final-state radi- 
ation. 



A. Partial Branching Fractions 

We correct the measured q^ spectra for resolution, ra- 
diative effects and bremsstrahlung by applying an un- 
folding technique that is based on singular- value decom- 
position of the detector response matrix [53. The de- 
tector response matrix in the form of a two-dimensional 
histogram of the reconstructed versus the true q^ values 
(see Figure 2]) is used as input to the unfolding algo- 
rithm. This algorithm contains a regularization term to 
suppress spurious oscillations originating from statistical 
fluctuations. To find the best choice of the regulariza- 
tion parameter k, we have studied the systematic bias on 
the partial branching fractions compared to the statisti- 
cal uncertainty as a function of k, using a set of simulated 
distributions. The data samples in this analysis are large 
enough that no severe distortions due to statistical fluc- 
tuations are expected. We choose the largest possible 
value of K, i.e., we set k equal to the number of q^ bins, 
to minimize a potential bias. 

The AS/Ag^ distributions resulting from the unfold- 
ing procedure are presented in Figure[21]for B -^ nlv and 
in Figure [25] for B ^ piu. Tables 1x1 and IXj list the par- 
tial branching fractions AS for B — >■ iilv and B — >■ ptv, 
respectively. 



B. Form- factor Shape 

For B — > Tilv decays, we extract the shape of the form 
factor f+{q^) directly from data. For B -^ piv decays, 
we restrict ourselves to the measurement of the q^ de- 
pendence, since the current experimental precision is not 
adequate to extract the three different form factors in- 
volved. 

Several parameterizations of f+{q^) are used to inter- 
polate between results of various form-factor calculations 
or to extrapolate these calculations from a partial to the 
whole q^ range. The four most common parameteriza- 
tions, the BK [H, BZ [H, BGL [H, HI and BCL 
parameterizations, have been introduced in Section [ill 
For the BGL and BCL parameterizations, we consider a 
linear {kmax = 2) and a quadratic {kmax = 3) ansatz. 



We perform x^ fits to the measured q^ spectrum to 
determine the free parameters for each of these parame- 
terizations. The fit employs the following x^ definition, 
with integration of the fit function over the q^ bins. 



Wbin. 



E ^^ ^7'^^ ' 



(32) 



i.j=i 



where V^ ■ is the inverse covariance matrix of the partial- 
branching- fraction measurements. A^ for q- bin k is de- 
fined as 



{AB\ '^^"' C 



VAgV, Ag2 



j pl\U{q^;a)\Hq\ (33) 



^il 



where a denotes the set of parameters for a chosen pa- 
rameterization of !+{q^), and C ~ |Kt{,proG|./(247r^) is 
an overall normalization factor whose value is irrelevant 
for these fits since the data can only constrain the shape 
of the form factor, but not its normalization. 

In Table IXIII and Figure [24] we present the results of 
these fits to the B — >■ irli/ samples. All parameterizations 
describe the data well, with x^ probabilities ranging from 
10% to 18%. Thus, within the current experimental pre- 
cision, all parameterizations are valid choices, and the 
central values for |Ktb|/+(0) agree with each other. We 
choose the quadratic BGL parameterization as the de- 
fault, though even a linear parameterization results in a 
very good fit to the data. The error band represents the 
uncertainties of the fit to data, based on the quadratic 
BGL parameterization (solid line in Figure [M)) . It has 
been computed using standard error propagation, taking 
the correlation between the fit parameters into account. 

We compare the measured q^ spectra with the shapes 
predicted by form-factor calculations based on lat- 
tice QC p [23j . light-cone sum rules [1^ [I^l, and the 
ISGW2 [ij] relativistic quark model. Among the avail- 
able calculations for B — > ttIi/ decays, the HPQCD lat- 
tice calculation agrees best with the data. It should 
be noted that the LQCD predictions are only valid for 
q^ > 16GeV^ the earlier LCSR calculation (LCSR 1) 
for g^ < 16GeV^, and the more recent LCSR calcula- 
tion (LCSR 2) for q'^ < 12GeV^; their extrapolation is 
impacted by sizable uncertainties. 

In Table IXll and Figure [25| we present the results of the 
fits to the B — > pii/ samples. The LCSR calculation and 
the ISGW2 model are in good agreement with the data. 
However, the errors of the measured B — > plv partial 
branching fractions are relatively large, at the level of 
15-30%, depending on the q"^ interval. 

It should be noted that the theoretical calculations 
differ most for low and high q^. In these regions of 
phase space, the measurements are impacted significantly 
by higher levels of backgrounds, specifically continuum 
events at low q^ and other B — >■ X^tv decays that are 
difficult to separate from the signal modes at higher 
q^ . These two background sources have been examined 
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TABLE X: Partial and total branching fractions (corrected for radiative effects) for B^' — !> Tv^l'^h' and B^ — )> tt^^+j/ decays 
obtained from the single-mode fits and B — >■ ttIv decays from the four-mode fit with statistical (fit), systematic and total errors. 
The branching fraction for B^ — ^ n^l'^u has been scaled by twice the lifetime ratio of neutral and charged B mesons. All 
branching fractions and associated errors are given in units of 10~*. 



q^ range (GeV^) 


0-4 


4-8 


8-12 


12-16 


16-20 


20-26.4 


Total 


<16 


>16 


AB{B° -^ -R-l+u) 


0.313 


0.329 


0.241 


0.222 


0.206 


0.124 


1.435 


1.105 


0.330 


Fit error 


0.030 


0.018 


0.018 


0.020 


0.020 


0.018 


0.061 


0.049 


0.027 


Syst. error 


0.025 


0.016 


0.015 


0.015 


0.013 


0.010 


0.068 


0.059 


0.019 


Total error 


0.039 


0.024 


0.023 


0.025 


0.024 


0.021 


0.092 


0.077 


0.033 


AS(B+ -^ -k"£+u) X 2ro/r+ 


0.357 


0.294 


0.234 


0.210 


0.206 


0.099 


1.401 


1.096 


0.305 


Fit error 


0.049 


0.031 


0.031 


0.033 


0.039 


0.043 


0.102 


0.075 


0.062 


Syst. error 


0.050 


0.015 


0.028 


0.019 


0.024 


0.028 


0.106 


0.089 


0.037 


Total error 


0.070 


0.035 


0.041 


0.038 


0.046 


0.051 


0.147 


0.117 


0.072 


AB{B° -^ -R-l+u) 4-mode 


0.320 


0.321 


0.235 


0.220 


0.201 


0.118 


1.414 


1.095 


0.319 


Fit error 


0.025 


0.017 0.015 


0.017 


0.018 


0.016 


0.050 


0.041 


0.024 


Syst. error 


0.027 


0.012 


0.016 


0.018 


0.014 


0.014 


0.074 


0.061 


0.024 


Total error 


0.037 


0.021 


0.022 


0.025 


0.023 


0.022 


0.089 


0.074 


0.034 
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FIG. 24: Measured AS/Ag^ distribution for B — >■ -kIv. The vertical error bars correspond to the combined statistical and 
systematic uncertainties. The positions of the data points have been adjusted to correspond to the mean q^ value in each bin, 
based on the quadratic BGL ansatz. Left: fits of four different form-factor parameterizations to the AS/Ag data spectrum. 
The fit result for the BZ and BCL parameterizations are barely visible, since they overlap almost completely with the BGL 
result. The shaded band illustrates the uncertainty of the quadratic BGL fit to data. Right: shape comparisons of the data to 
various B -> -Rtv form-factor predictions (LCSR 1 [H, LCSR 2 [Til, HPQCD [H, ISGW2 [11]), which have been normalized 
to the measured total branching fraction. The extrapolations of the QCD predictions to the full q^ range are marked as dashed 
lines. 



in detail, and the uncertainties in their normalization 
and shape are included in the systematic uncertainties. 
For the inclusive B -^ X^^v background, the q^ and 
the hadronic mass spectra arc derived from theoretical 
predictions that depend on non-perturbative parameters 
that are not well measured [Sg. For B — > plv the corre- 
lation between the signal and the B — >■ X^tv background 
is so large that they cannot both be fitted simultaneously. 



Thus the B — ?> X^^v background scale factor and shape 
are fixed to the MC predictions, which have large uncer- 
tainties. MC studies indicate that this may introduce a 
bias affecting the signal yield. The stated errors account 
for this potential bias. 
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TABLE XI: Partial and total branching fractions (corrected 
for radiative effects) for B'^ — >■ p~ l^v and B^ — >■ fPt^v decays 
obtained from the single-mode fits and for B — >■ plv decays 
from the four- mode fit with statistical (fit), systematic and 
total errors. The branching fractions for B^ — >■ fPl^v have 
been scaled by twice the ratio of the lifetimes of neutral and 
charged B mesons. All branching fractions and associated 
errors are given in units of 10"*. 



q^ range (GeV^) 


0-8 8-16 


16-20.3 


Total 


AZ3(B" -^ p- 


-e+v) 




0.747 0.980 


0.256 


1.984 


Fit error 






0.151 0.087 


0.030 


0.214 


Syst. error 






0.178 0.165 


0.066 


0.379 


Total error 






0.234 0.187 


0.072 


0.435 


Ae(B+ ^ p 


H+v) 


X 2ro/r+ 


0.627 0.977 


0.265 


1.871 


Fit error 






0.136 0.079 


0.028 


0.190 


Syst. error 






0.152 0.161 


0.061 


0.320 


Total error 






0.204 0.179 


0.068 


0.373 


Ae(B° -^ p- 


-tv) 


4-mode 


0.564 0.912 


0.268 


1.745 


Fit error 






0.107 0.059 


0.022 


0.149 


Syst. error 






0.126 0.135 


0.058 


0.272 


Total error 






0.166 0.147 


0.062 


0.310 
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FIG. 25: (color online) Measured ^Bj t\c^ distribution for 
B — >■ piv. The inner and outer error bars correspond to 
the statistical uncertainty and the combined statistical and 
systematic uncertainty, respectively. The data are compared 
with the B — >■ piv form-factor predictions from LCSR [131 
and from the ISGW2 quark model [11. 



measured partial branching fractions with integrals of the 
form-faetor calculations over a certain q^ range using the 
relation 



IKbl = 



A6(g: 



2 






To ^C{qLn^q?nax) 



(34) 



where tq = (1.530 ± 0.009) ps is the B° lifetime and AC 
is defined as 



^ V^mm ' qmax) 



Lr p 

247r3 



Pi\f+{q')\'dq' . (35) 



The values of A^ are derived from theoretical form-factor 
calculations for different q^ ranges. Table IXIIII summa- 
rizes the AC values, the partial branching fractions and 
the \Vub\ results. 

For B — >■ p£v, values of AC are taken from the LCSR 
calculation in the range q^ < 16 GeV and the quark model 
predictions of ISGW2 over the full q^ range. The results 
are also presented in Table IXllIl Estimates of the uncer- 
tainties for AC are not given in Refs. [131 ^^^ H ■ 

Second, we perform a simultaneous fit to the most re- 
cent lattice results and BaBAR data to make best use of 
the available information on the form factor from data 
(shape) and theory (shape and normalization). A fit 
of this kind was first presented by the FNAL/MILC 
Collaboration [23l using the earlier BaBAR results on 
B° -^ TT-^+v decays Q. 

To perform this fit, we translate the f+{q^) predictions 
from LQCD to l/(To|K6)P)AS/Ag2. We simultaneously 
fit this distribution and the AB/Aq^ distribution from 
data as a function of q^. We use the BGL form- factor 
parameterization as the fit function, with the additional 
normalization parameter anorm = TolVubp, which allows 
us to determine |V^6| from the relative normalization of 
data and LQCD predictions. 

The x^ for this fit is given by 

X^ = x^ (data) + x^ (lattice) 

\ ^ A data /T^data\ — 1 A data 

ATpointB 

+ Y. ^J'' ivt'r'^t' (36) 

/,m— 1 



C. Determination of |V^b| 

We choose two different approaches to determine the 
magnitude of the CKM matrix element Vub- 

First, we use the traditional method to derive |14f,|- 
As in previous publications [y, l8l-KLl|. we combine the 



where 



I data 



L lat 






G 



data 



1 

A^ 



g{q^;a)dq^ , (37) 



Ag2 



F Ji r„2\\ j-lat/„2\|2 



247r3 



Pi{qt)\fr{qtr-9{qt-.a) 



(38) 



35 



TABLE XII: Results of fits to the measured AB/Aq^ for B — > tvIu decays, based on different form-factor parameterizations. 



Parametrization x /ndf Prob(x /ndf) Fit parameters 



/+(0)|Kbl [10- 



BK 


6.8/4 


0.148 


BZ 


6.0/3 


0.112 


BCL (2 par.) 


6.3/4 


0.179 


BCL (3 par.) 


6.0/3 


0.112 


BGL (2 par.) 


6.6/4 


0.156 


BGL (3 par.) 


6.3/3 


0.100 



asK = +0.310 ± 0.085 
TBZ = +0.170 ±0.124 
asz =+0.761+0.337 
61/60 = -0.67 + 0.18 
61/60 = -0.90 + 0.46 
62/60 = +0.47+1.49 
ffli/tso = -0.94 + 0.20 
ffli/ao = -0.82 + 0.29 
02/^0 = -1.14+1.81 



1.052 + 0.042 

1.079 + 0.046 

1.065 + 0.042 
1.086 + 0.055 

1.103 + 0.042 

1.080 + 0.056 



TABLE XIII: \Vub\ derived from B -)■ ■kIu and B -^ plv 
decays for various q^ regions and form- factor calculations. 
Quoted errors are experimental uncertainties and theoretical 
uncertainties of the form-factor integral A^". No uncertainties 
on AC, for B — >■ plv are given in Refs. [131 and [1J|- 





q^ Range 


AB AC 


IKbl 




(GeV^) 


(10-*) (ps-i) 


(10-3) 


B-> Tilv 








LCSR 1 [15] 


0-16 


1.10 ± 0.07 5.44+1.43 


3.63 ± 0.12tgjg 


LCSR 2 [19] 


0-12 


0.88 + 0.06 4.00lJ;°^ 


3.78 + 0.13to:4o 


HPQCD [2^ 


16 - 26.4 


0.32 + 0.03 2.02+0.55 


3.21 + 0.17t;;;f6 


B-> plv 








LCSR [17] 


0-16.0 


1.48 + 0.28 13.79 


2.75 + 0.24 


ISGW2 [14] 


0-20.3 


1.75 + 0.31 14.20 


2.83 + 0.24 



and 



3(9^; a) 



Ui^) 



Onorm for data 



(39) 



1 for LQCD ' 



1 



r{q')<i^{q\ql) 



fe=0 



Here {AB / Aq'^)^'^^^ is the measured spectrum, f^^{qf) 
are the form-factor predictions from LQCD, and 
^ydntayi ^^^ (yUf-^-i g^j,g ^j^g corresponding in- 
verse covariance matrices for (AB/Aq^)*^^*^ and 
Gl/{247r^)pl{qf)\,ff{qf)\^, respectively. The set of free 
parameters a of the fit function g{q^; a) contains the co- 
efficients flfc of the BGL parameterization and the nor- 
malization parameter flnorm- 

From the FNAL/MILC ^ lattice calculations, we use 
only subsets with six, four or three of the twelve predic- 
tions at different values of q^, since neighboring points 
are very strongly correlated. All chosen subsets of LQCD 
points contain the point at lowest q^ . It has been checked 
that alternative choices of subsets give compatible re- 



sults. From the HPQCD [23| lattice calculations, we use 
only the point at lowest q^ since the correlation matrix 
for the four predicted points is not available. For com- 
parison, we also perform the corresponding fit using only 
the point at lowest q^ from FNAL/MILC. The data, the 
lattice predictions, and the fitted functions are shown in 
Figure [251 Table IXIVI shows the numerical results of the 
fit. 

For the nominal fit we use the subset with four 
FNAL/MILC points and assume a quadratic BGL pa- 
rameterization. We refer to this fit as 3+1-parameter 
BGL fit (three coefficients Ofc and the normalization pa- 
rameter anorm)- As Can be seen in Table [Xlll for the fit to 
data alone, the data are well described by a linear func- 
tion with the normalization ag and a slope ai/ao- This 
indicates that most of the variation of the form factor 
is due to well-understood QCD effects that arc parame- 
terized by the functions V{q^) and <f>{q'^,ql) in the BGL 
parameterization. If we include a curvature term in the 
fit, the slope ai/oo = —0.82 + 0.29 is fully consistent with 
the linear fit; the curvature aija^ is negative and consis- 
tent with zero. Since the z distribution is almost linear, 
we also perform a linear fit (2+1-parameter BGL fit) for 
comparison. The results of the linear fits are also shown 
in Table [XlVl 

The simultaneous fits provide very similar results, both 
for the BGL expansion coefficients, which determine the 
shape of the spectrum, and for jV^hj. The fitted values for 
the form-factor parameters are very similar to those ob- 
tained from the fits to data alone. This is not surprising, 
since the data dominate the fit results. Unfortunately 
the decay rate is lowest and the experimental errors are 
largest at large q^ , where the lattice calculation can make 
predictions. We obtain from these simultaneous fits 



\Yuh\ = (2.87+0.28) X IQ-^ 

\Yuh\ ^ (2.95 + 0.31) X 10-3 

\Vub\ = (2.93 + 0.31) X 10-3 

\Vub\ = (2.92 + 0.37) X 10-3 

\Vub\ = (2.99 + 0.35) X 10-3 



FNAL/MILC (6 points) 
FNAL/MILC (4 points) 
FNAL/MILC (3 points) 
FNAL/MILC (1 point) , 
HPQCD (1 point) , 



where the stated error is the combined experimental and 
theoretical error obtained from the fit. The coefficients 
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FIG. 26: (color online) Simultaneous fits of the BGL parameterization to data (solid points with vertical error bars representing 
the total experimental uncertainties) and to four of the twelve points of the FNAL/MILC lattice prediction (magenta, closed 
triangles). Left: linear (2+1-parameter) BGL fit, right: quadratic (3+1-parameter) BGL fit. The LQCD results are reseated 
to the data according to the \Vub\ value obtained in the fit. The shaded band illustrates the uncertainty of the fitted function. 
For comparison, the HPQCD (blue, open squares) lattice results are also shown. They are used in an alternate fit. 



Ofc are significantly smaller than 1, as predicted. The 
sum of the squares of the first two coefhcients, J2k=o "fc ~ 
(0.85±0.20) X 10""^, is consistent with the tighter bounds 
set by Becher and Hill [25i] . 

Since the total error of 10% on \Vub\ results from the 
simultaneous fit to data and LQCD predictions, it is non- 
trivial to separate the error into contributions from ex- 
periment and theory. We have estimated that the error 
contains contributions of 3% from the branching-fraction 
measurement, 5% from the shape of the q^ spectrum de- 
termined from data, and 8.5% from the form-factor nor- 
malization obtained from theory. 

We study the effect of variations of the isospin rela- 
tions imposed in the combined four-mode fit as stated 
in Eqs. [301 These relations are not expected to be ex- 
act, though the comparison of the single-mode fit results 
provides no indication for isospin breaking. The isospin- 
breaking effects are primarily due to tt" — t; and p" — w 
mixing in B+ — > 7r°£+zy and 5+ -^ p^i'^v decays, re- 
spectively. They are expected to increase the branching 
fractions of the B'^ relative to the B^ meson. Given the 
masses and widths of the mesons involved, the impact 
of TT*^ — ?7 mixing is expected to be smaller than that of 
p" — a; mixing. 

Detailed calculations have been performed to cor- 
rect form-factor measurements and to extract Vus from 
semileptonic decays of charged and neutral kaons [54| . 
These calculations account for isospin breaking due to 
tt" — ?7 mixing and should also be applicable to B~^ — >■ 



7r°^+zy decays. For B^ 



rOp+, 



decays the effect is 



expected to be smaller by a factor of three, i.e., the pre- 
dicted increase is (1.5 ± 0.2)% [5|. For B+ -^ p°(.+v de- 
cays, calculations have not been carried out to the same 
precision. Based on the change in the tt+tt" rate at the 



peak of the p mass distribution, the branching fraction 
is predicted to increase by as much as 34% [56[- How- 
ever, an integration over the resonances weighted by the 
proper Breit-Wigner function and taking into account 
the masses and finite p and to widths results in a much 
smaller effect, an increase in the tt+tt^ branching fraction 

of 6% i^. 



We have assessed the impact of changes in the ratios of 
the branching fractions for charged and neutral B mesons 
on the extraction of the differential decay rates due to ad- 
justments of the MC default branching fractions of the 
5+ decays in the combined four-mode fit. For a 1.5% in- 
crease in the i3+ — >■ ir^t.'^v branching fraction, the fitted 
B — > niv partial branching fraction decreases by 0.5%, 
while the B — >■ p(.v rate increases by less than 0.1%. A 
6% increase in the B~^ — > p^i'^v branching fraction re- 
sults in a decrease of the B — >■ piv rate by 3.1% and a 
0.14% increase for the fitted B — )• irlv rate. We observe 
a partial compensation to the change in the simulated 



5+ — > -K^t^v rate due to changes in the B^ 



P 



H+u 



background contribution, and vice versa. The observed 
changes in the fitted yields depend linearly on the im- 
posed branching-fraction changes and are independent 
of q^. 



For a 1.5% variation of the i?+ — > n^i'^v branching 
fraction, the value for \Vub\ extracted from the measured 
B -^ ■niv spectrum decreases by 0.2%. A -1-6% variation 
of the B^ -^ p^l^v branching fraction increases the value 
of \Vub\ extracted from the same measured spectrum by 
0.3%. 
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TABLE XIV: Results of simultaneous fits to data and LQCD calculations, based on the linear or quadratic BGL parameteri- 
zations. \Vub\ is determined from the relative normalization of data and theory prediction. 



Parameterization 




Input 


xVndf Prob(xVndf) 


Fit parameters 


BGL (2+1 par.) 




Data 


18.2/9 


0.033 


ao = (2.07 + 0.21) x 10"^ 




+ 


FNAL/MILC 
(6 points) 






ai/ao = -0.78 + 0.22 

\Vub\ = (3.04 + 0.38) X 10-^ 


BGL (2+1 par.) 




Data 


7.1/7 


0.415 


ao = (2.16 + 0.19) X lO-"* 




+ 


FNAL/MILC 
(4 points) 






ai/ao = -0.93 + 0.20 

\Vub\ = (2.99 + 0.32) X 10'^ 


BGL (2+1 par.) 




Data 


6.8/6 


0.341 


ao = (2.17 + 0.19) X 10"^ 




+ 


FNAL/MILC 
(3 points) 






ai/ao = -0.93 + 0.20 

jKfcl = (2.97 + 0.32) X 10"^ 


BGL (2+1 par.) 




Data 


6.6/4 


0.156 


ao = (2.23 + 0.26) x 10"^ 




+ 


FNAL/MILC 
(1 point) 






ai/ao = -0.94 + 0.20 

\Vub\ = (2.90 + 0.36) X 10-^ 


BGL (2+1 par.) 




Data 

+ HPQCD 

(1 point) 


6.6/4 


0.156 


ao = (2.19 + 0.23) x lO^^ 

ai/ao = -0.94 + 0.20 

\Vub\ = (2.94 + 0.34) X 10"^ 



BGL (3+1 par.; 



BGL (3+1 par.) 



BGL (3+1 par.; 



BGL (3+1 par. 



BGL (3+1 par.; 



Data 9.8/8 

FNAL/MILC 
(6 points) 



0.276 



ao = (2.31 + 0.20) x 10"^ 
ai/ao = -0.71 + 0.19 
a2/ao = -2.33 + 0.84 



IVub 



(2.87 + 0.28) X 10" 



Data 6.6/6 

FNAL/MILC 
(4 points) 



0.355 ao = (2.22 + 0.21) x 10"^ 

ai/ao = -0.86 + 0.23 
a2/ao = -0.97+1.36 
IKtl = (2.95 + 0.31) X 10"^ 



Data 6.3/5 

FNAL/MILC 
(3 points) 



0.279 ao = (2.24 + 0.22) x 10"^ 

ai/ao = -0.84 + 0.23 
a2/ao = -1.01 + 1.40 
\Vub\ = (2.93 + 0.31) X 10"^ 



Data 6.3/3 

FNAL/MILC 
(1 point) 



0.100 ao = (2.24 + 0.26) x 10"^ 

ai/ao = -0.82 + 0.29 
a2/ao = -1.14 + 1.81 
\Vub\ = (2.92 + 0.37) X 10"^ 



Data 
+ HPQCD 

(1 point) 



6.3/3 0.100 ao = (2.19 + 0.23) X 10"^ 

ai/ao = -0.82 + 0.29 
a2/ao = -1.14 + 1.81 
\Vub\ = (2.99 + 0.35) X 10"^ 



IX. CONCLUSIONS 

In summary, we have m.easured the exclusive branching 
fractions B{B^ -^ tt-(.+v) and B{B° -^ p-(.+v) as a 
function of q^ and have determined \Vub\ using recent 
form-factor calculations. We measure the total branching 
fractions, based on samples of charged and neutral B 
mesons and isospin constraints, to be 



6(5" 
6(B" 



■ p-t 



-y) = 



= (1.41 ±0.05 ±0.07) X lO"'' 
>) = (1.75 ±0.15 ±0.27) X 10 



-4 



where the first error is the statistical uncertainty of the fit 
employed to determine the signal and background yields 



and the second is the systematic uncertainty. The sepa- 
rate measurements of the branching fractions for charged 
and neutral B mesons arc consistent within errors with 
the assumed isospin relations. 



g(B" ^ 7r^£+i/) 
B{B+ -^ TT^i+v) X 2^ 

BjB^ -^ p-l+v) 
B{B+ -^ p°£+iy) X 2^ 



1.03 ±0.09 ±0.06 , 
1.06 ±0.16 ±0.08 . 



We have assessed the sensitivity of the combined 
branching-fraction measurements to isospin violations 
due to TT^ — r] and p^ — to mixing in B^ decays. Based 
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on the best estimates currently available, the impact on 
the branching fractions is small compared to the total 
systematic errors. We refrain from applying corrections, 
given the uncertainties in the size of the effects. 

The measured branching fraction for B — > iriv is 
more precise than any previous measurement and agrees 
well with the current world average B{B^ — >■ 7r~£+j/) = 
(1.36 ± 0.05 ± 0.05) X 10-4 0]. The branching fraction 
for B — >■ piv is also the most precise single measurement 
to date based on a large signal event sample, although 
the Belle Collaboration [il| has reported a smaller sys- 
tematic error (by a factor of two), based on a small signal 
sample of hadronically-tagged events [11| . The B — >■ pHv 
branching fraction presented here is significantly lower 
(by about 2.5 a) compared to the current world aver- 



41 



age B{B° -^ p-i+u) = (2.77 ± 0.18 ± 0.16) x 10 
The dominant uncertainty of this B — > ptv measure- 
ment is due to the limited knowledge of the normaliza- 
tion and shape of the irreducible background from other 
B -^ Xu^v decays. 

Within the sizable errors, the measured q^ spectrum 
for B — >■ p^v agrees well with the predictions from light- 
cone sum rules [141 and the ISGW2 [1J| quark model. 
Neither of these calculations includes an estimate of their 
uncertainties. In the future, it will require much cleaner 
data samples and considerably better understanding of 
other B — >■ X^^-v decays to achieve significant improve- 
ments in the measurements of the form factors in B de- 
cays to vector mesons. 

For B — )- irtv decays, the measured q^ spectrum agrees 
best with the one predicted by the HPQCD lattice cal- 
culations [2j] . The measurement of the differential deca; 
rates is consistent with earlier BaBAR measurements [6|, 
within the stated errors, though the yield at low q^ is 
somewhat higher than previously measured. This results 
in a smaller value of as k, the parameter introduced by 
Becirevic and Kaidalov l2al, namely asK = 0.31 ± 0.09. 
Using the BGL ansatz, we determine a value /+(0)|Ku6| = 
(1.08 ± 0.06) X 10"^, which is larger than the value, 
/+ (0)1 Kb I == (0.91±0.06±0.3) X 10-3 ra, based on the 
earlier BaBAR decay rate measurement [9| and an average 
branching fraction of (1.37 ± 0.06 ± 0.07) x 10"'' [4l|. 

We determine the CKM matrix element \Vub\ using 
two different approaches. First, we use the traditional 
method to derive \Vub\ by combining the measured par- 
tial branching fractions with the form-factor predictions 
based on different QCD calculations. The results, pre- 
sented in Table IXIIIi agree within the sizable uncertain- 
ties of the form-factor predictions. For this approach we 
quote as a result the value of 

|Kb| = (3.78±0.13t^■j^)xl0-^ 



very well. We quote as a result the fitted value of 

iKbH (2.95 ± 0.31) X 10-3, 

based on the normalization predicted by the 
FNAL/MILC Collaboration. The total error of 

10% is dominated by the theory error of 8.5%. This 
value of \Vub\ is smaller by one standard deviation 
compared to the results of a combined fit to earlier 
BaBAR measurements and the same recent FNAL/MILC 
lattice calculations [23 |. 

The values of \Vub\ presented here appear to be sensi- 
tive to the q^ range for which theory predictions and the 
measured spectrum can be compared. LCSR calculations 
are restricted to low values of q^ and result in values of 
\Vub\ in the range of (3.63 — 3.78) x IQ-^ with theoretical 
uncertainties of _j^]^% and experimental errors of 3 — 4%. 
LQCD predictions are available for q^ > 16 GeV and 
result in |T4b| in the range of (2.95 — 3.21) x IQ-^ and 
experimental errors of 5 — 6% for both the traditional 
method and the simultaneous fit to LQCD predictions 
and the measured spectrum. This fit combines the mea- 
sured shape of the spectrum over the full q^ range with 
the lattice QCD form-factor predictions at high q^ and 
results in a reduced theoretical uncertainty of 8.5%, as 
compared to l}i% for the traditional method. 

Both \Vu,b\ values quoted as results are also lower 
than most determinations of \Vub\ based on inclusive 
B — > Xutv decays, which arc typically in the range 
(4.0 — 4.5) X 10-3. These inclusive measurements are 
very sensitive to the mass of the b quark, which is ex- 
tracted from fits to moments of inclusive B — > Xci-v and 
B ^ Xs"f decay distributions [33| and depends on higher- 
order QCD corrections. Estimated theoretical uncertain- 
ties are typically 6%. 

Global fits constraining the parameters of the CKM 
unitarity triangle performed by the CKMfitter [S^ and 
UTfit [591 Collaborations currently predict values for 
\Vub\ that fall between the two results presented here, 
\Vub\ = 3.5ltg;}^ X 10-3 and |Kh| = 3.41 ± 0.18 x 10-3, 
respectively. 

To permit more stringent tests of the CKM frame- 
work and its consistency with the standard model of 
elcctroweak interactions, further reductions in the ex- 
perimental and theoretical uncertainties will be neces- 
sary. For B — > ttIv decays this will require a reduction 
in the statistical errors and improved detector hermetic- 
ity to more effectively reconstruct the neutrino, which 
will reduce backgrounds from all sources. Further im- 
provements in the precision of lattice and other QCD 
calculations will also be beneficial. 



< 



12 GeV^. Second, we extract \Vub\ from simultaneous fits 
to data and lattice predictions using the quadratic BGL 
parameterization for the whole q'^ range. These fits to 
data and the two most recent lattice calculations by the 
FNAL/MILC [22I and HPQCD [H Collaborations agree 
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TABLE XV: Systematic errors in % for B{B" -^ n'l+f) from the B° -H- n'£+v and B+ -^ n^l+f one-mode fits. 









B 


^ ^ n-£+,. 








B 


+ ^ n^£+.. 






q^ range (GeV^) 


0-4 4-8 


8-12 


12-16 


16-20 


>20 0-26.4 


0-4 


4-8 


8-12 


12-16 


16-20 


>20 0-26.4 


Track efficiency 


2.0 


1.7 


2.9 


1.3 


0.1 


2.3 


1.7 


7.3 


1.8 


3.3 


1.6 


1.8 


6.8 


3.7 


Photon efficiency 


0.7 


0.9 


1.9 


3.5 


0.7 


1.0 


1.6 


2.9 


1.8 


5.3 


2.4 


8.9 


10.2 


4.8 


Lepton identification 


4.1 


1.7 


1.8 


1.8 


2.1 


2.4 


1.8 


3.3 


1.3 


2.2 


1.4 


2.7 


2.8 


1.4 


Kl efficiency 


1.0 


0.1 


0.3 


3.5 


1.2 


1.7 


1.2 


1.2 


0.5 


2.4 


3.6 


1.8 


1.1 


1.8 


Kl shower energy 


0.1 


0.1 


0.1 


0.7 


1.4 


2.8 


0.7 


2.6 


0.5 


0.4 


0.6 


1.4 


4.1 


1.4 


Kl spectrum 


1.1 


2.0 


2.5 


2.6 


5.1 


1.1 


2.4 


2.7 


0.7 


2.0 


4.5 


4.4 


4.2 


2.9 


B -^ ■Kf.vFF f+ 


0.5 


0.3 


0.4 


0.5 


0.6 


0.8 


0.5 


0.5 


0.5 


0.6 


0.7 


1.8 


2.6 


1.0 


B -^ pluFFAi 


1.6 


1.8 


2.0 


0.9 


2.0 


1.5 


1.6 


2.7 


2.2 


3.6 


2.0 


1.6 


5.1 


2.8 


B -i> pivFFA2 


1.5 


1.7 


1.7 


0.5 


2.1 


0.6 


1.4 


2.5 


2.0 


3.0 


1.4 


1.9 


2.7 


2.2 


B -^ plvFFV 


0.3 


0.4 


0.5 


0.3 


0.4 


0.8 


0.4 


1.1 


0.7 


1.2 


1.0 


0.9 


3.0 


1.2 


B{B" -^ p-l+v) 


0.3 


0.2 


0.5 


0.8 


0.7 


0.4 


0.5 


0.4 


0.2 


0.3 


0.4 


0.9 


1.6 


0.5 


B(B+ -^ Lue+iy) 


0.1 


0.1 


0.1 


0.1 


0.3 


1.4 


0.2 


0.6 


0.3 


0.6 


0.6 


0.9 


4.8 


1.1 


B{B+ -^ r)l+v) 


0.1 


0.1 


0.1 


0.1 


0.1 


0.1 


0.1 


0.1 


0.1 


0.1 


0.1 


0.5 


0.6 


0.2 


B{B+ -> r)'l+v) 


0.1 


0.1 


0.1 


0.1 


0.1 


0.3 


0.1 


0.1 


0.1 


0.1 


0.1 


0.4 


0.7 


0.2 


B{B ^ Xuli^) 


0.1 


0.1 


0.1 


0.1 


0.3 


0.9 


0.2 


0.7 


0.2 


0.3 


0.4 


1.2 


1.4 


0.6 


B -^ Xu.lv SF param 


0.4 


0.2 


0.1 


0.4 


0.6 


2.3 


0.5 


0.6 


0.3 


0.4 


1.1 


1.3 


6.2 


1.3 


B -> Dlv FF p^ 


0.1 


0.1 


0.4 


0.3 


0.1 


0.3 


0.2 


0.4 


0.1 


0.6 


0.3 


0.2 


0.8 


0.4 


B -^ D'lu FF Ri 


0.2 


0.4 


0.8 


0.9 


0.3 


0.9 


0.6 


0.1 


0.5 


1.3 


0.5 


0.7 


1.8 


0.7 


B -^ D'tv FF R2 


0.5 


0.1 


0.2 


0.1 


0.1 


0.2 


0.2 


0.5 


0.3 


0.6 


0.5 


0.9 


1.4 


0.6 


B -^ D'lu FF p]^. 


0.9 


0.4 


0.5 


1.0 


0.2 


0.8 


0.6 


0.7 


0.4 


0.6 


0.6 


0.6 


2.4 


0.8 


B{B -^ D£iy) 


0.2 


0.1 


0.2 


0.3 


0.4 


0.1 


0.2 


1.1 


0.3 


0.8 


0.3 


0.3 


0.9 


0.6 


B{B -^ D*lv) 


0.5 


0.1 


0.3 


0.4 


0.3 


0.3 


0.3 


0.9 


0.4 


0.6 


0.7 


1.7 


1.8 


0.9 


B{B -^ i3**^l/)narrow 


0.3 


0.1 


0.1 


0.4 


0.1 


0.1 


0.2 


0.2 


0.2 


0.3 


0.3 


0.2 


0.2 


0.2 


B{B ^ D"£i/)taoad 


0.2 


0.1 


0.1 


0.5 


0.1 


0.2 


0.2 


0.3 


0.3 


0.3 


0.4 


0.1 


0.6 


0.3 


Secondary leptons 


0.8 


0.1 


0.2 


0.1 


0.1 


0.2 


0.2 


0.5 


0.2 


0.3 


0.4 


0.8 


0.9 


0.5 


Continuum 


5.0 


1.6 


1.6 


1.3 


1.8 


4.1 


1.7 


8.8 


1.8 


6.1 


4.0 


4.8 


13.1 


4.2 


Bremsstrahlung 


0.2 


0.1 


0.1 


0.1 


0.1 


0.2 


0.1 


0.2 








0.1 


0.2 


0.2 


0.1 


Radiative corrections 


0.4 


0.1 


0.1 


0.3 


0.5 


0.9 


0.3 


0.5 


0.1 


0.1 


0.3 


0.5 


0.9 


0.4 


Nbb 


1.1 


1.2 


1.1 


1.2 


1.2 


1.0 


1.1 


1.2 


1.1 


1.2 


1.2 


1.1 


1.2 


1.2 


f±/foo 


1.5 


1.3 


1.3 


1.3 


1.0 


0.8 


1.3 


1.3 


1.1 


1.0 


1.0 


0.9 


0.2 


1.0 


Total 


7.8 


4.8 


6.0 


6.8 


7.1 


7.6 


4.4 


13.7 


5.0 


11.1 


8.6 


12.7 


22.4 


7.5 
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TABLE XVI: Systematic errors in % for B{B° -)> p'l+i^) from tlie B° -5> p~£+iy and B+ -> p"^+i^ one-mode fits. 





B^^ 


p-l+v 


B+ ^ p^e+i. 


g" range (GeV^) 


0-8 


8-16 


> 16 0-20.3 


0-8 


8-16 


> 16 0-20.3 


Track efficiency 


2.4 


0.8 


3.0 


1.8 


2.4 


4.0 


1.0 


2.9 


Photon efficiency 


1.4 


1.5 


3.9 


2.0 


5.5 


3.5 


3.5 


4.2 


Lepton Identification 


2.6 


2.9 


4.7 


3.0 


4.3 


3.2 


4.6 


3.7 


Kl efficiency 


8.2 


1.0 


6.5 


4.3 


10.1 


1.8 


5.8 


5.2 


Kl shower energy 


4.4 


0.2 


0.4 


1.6 


1.0 


0.7 


0.8 


0.8 


Kl spectrum 


9.3 


8.2 


9.7 


8.8 


1.0 


5.2 


7.3 


4.3 


B -^ Tilv FF /+ 


1.2 


0.3 


2.8 


1.1 


0.8 


0.2 


2.4 


0.8 


B -5> plv FF yli 


14.8 


8.3 


4.8 


9.6 


14.3 


8.3 


4.3 


9.4 


B -^ p£u FF A2 


11.3 


5.1 


0.8 


6.2 


11.0 


5.0 


0.6 


6.0 


B -^ plf FF V 


3.7 


3.5 


3.9 


3.6 


3.8 


3.6 


3.5 


3.6 


B{B" -^ n-£+iy) 


0.8 


0.6 


1.8 


0.9 


0.6 


0.1 


1.5 


0.5 


B{B+ -f Lol+v) 


0.1 


0.4 


3.3 


0.8 


0.9 


0.9 


3.3 


1.4 


B{B+ -^ -ql+v) 


0.1 


0.1 


0.8 


0.2 


0.1 


0.1 


0.6 


0.1 


B{B+ -^ rj'£+u) 


0.1 


0.1 


1.2 


0.3 


1.2 


0.8 


1.1 


1.0 


B{B -^ Xu£u) 


4.0 


5.3 


10.3 


5.9 


7.9 


8.0 


11.0 


8.6 


B — ;> Xu£v SF param. 


6.9 


6.5 


13.1 


7.9 


10.9 


7.8 


12.7 


9.7 


B -^ D£v FF pi:, 


0.1 


0.2 


0.1 


0.1 


0.6 


0.1 


0.1 


0.2 


B -5> D''£v¥¥ Ri 


0.1 


0.1 


0.2 


0.1 


0.1 


0.1 


0.1 


0.1 


B -> D*£uFF R2 


0.2 


0.1 


0.3 


0.2 


0.6 


0.1 


0.2 


0.3 


B -^ D*euFF pI. 


0.4 


0.2 


0.1 


0.2 


1.4 


0.1 


0.2 


0.5 


B{B -^ D£v) 


0.6 


0.3 


0.2 


0.4 


0.4 


0.2 


0.2 


0.3 


B{B -^ D*£iy) 


0.3 


0.1 


0.2 


0.2 


0.2 


0.1 


0.2 


0.2 


B{B -^ D**£u)n.rro^ 


0.7 


0.1 


0.1 


0.2 


0.6 


0.1 


0.1 


0.2 


B{B -)■ D**^!/)broad 


0.9 


0.1 


0.1 


0.3 


0.9 


0.1 


0.1 


0.4 


Secondary leptons 


0.5 


0.1 


0.2 


0.2 


0.7 


0.1 


0.1 


0.2 


Continuum 


5.9 


3.4 


6.3 


3.2 


8.2 


3.8 


6.4 


3.4 


Bremsstrahlung 


0.2 


0.1 


0.2 


0.2 


0.2 


0.1 


0.2 


0.2 


Radiative corrections 


0.6 


0.1 


0.8 


0.4 


0.6 


0.1 


0.8 


0.4 


Nbb 


1.7 


1.9 


2.9 


2.0 


2.2 


2.1 


2.9 


2.3 


/±//oo 


1.9 


1.4 


1.0 


1.5 


1.6 


1.4 


0.7 


1.3 


Total 


25.6 


16.5 


23.9 


19.4 


27.6 


18.0 


22.5 


17.1 
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B. Correlation and Covariance Matrices 

Table IXVIII shows the fuU correlation matrix for all 
signal and background fit parameters in the four-mode 
maximum-likelihood fit used to determine the signal 
yields, described in Section |Vll This appendix also con- 
tains all statistical, systematic and total correlation and 
covariance matrices for the B — ^ niiy and B — > piv 
AB/Aq^ measurements. The total correlation matrix is 
shown before and after unfolding of the q^ spectrum. All 
covariance matrices are shown after q^ unfolding. The 
total covariance matrix for B — > ttEv in Table IXXVIIll is 
used in the form- factor fits described in Eq. |31] or [351 





1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


14 


15 


16 


17 


18 


19 


20 


21 


22 


23 


24 25 


pV± 1 


1.000 
















































p^±' 


' 2 


0.147 


1.000 














































otherBB o 


-0.654 


-0.449 


1.000 












































p:±, 


i4 


-0.096 


-0.449 


-0.225 


1.000 










































p:±. 


.5 


-0.175 


-0.119 


0.033 


0.214 


1.000 








































pf 


' 6 


0.021 


0.197 


-0.148 


-0.158 


-0.080 


1.000 






































pf 


■" 7 


0.055 


0.060 


-0.066 


-0.158 


-0.111 


0.264 


1.000 




































pf 


' 8 


0.033 


0.002 


-0.015 


-0.079 


-0.564 


0.137 


0.189 


1.000 


































pT* 


' 9 


-0.263 


0.134 


0.030 


-0.013 


0.044 


-0.144 


-0.033 


-0.019 


1.000 
































P^" 


10 


-0.105 


0.126 


-0.007 


-0.120 


0.032 


-0.429 


-0.050 


-0.033 


0.191 


1.000 






























Pj^" 


11 


0.086 


0.150 


-0.110 


-0.271 


-0.044 


0.085 


-0.156 


0.024 


0.050 


0.089 


1.000 




























Pf" 


12 


0.004 


-0.073 


0.110 


-0.232 


-0.025 


-0.001 


-0.267 


0.026 


-0.005 


0.058 


0.197 


1.000 


























vf- 


13 


-0.162 


0.116 


0.244 


-0.432 


0.216 


0.061 


0.059 


-0.469 


0.068 


0.085 


0.127 


0.135 


1.000 
























Pi"- 


14 


-0.124 


0.036 


0.079 


-0.033 


-0.602 


0.007 


-0.005 


0.035 


0.057 


0.011 


0.005 


-0.008 


0.032 


1.000 






















P% 


15 


0.055 


-0.015 


-0.049 


0.029 


0.103 


0.031 


0.037 


-0.020 


-0.119 


0.059 


-0.050 


0.056 


-0.079 


-0.163 


1.000 




















p"J'" 


16 


-0.075 


0.104 


-0.009 


-0.082 


-0.043 


0.128 


-0.019 


0.024 


0.156 


0.109 


0.086 


-0.021 


0.095 


0.073 


-0.193 


1.000 


















^..H.-rBB 


17 


0.003 


-0.086 


0.092 


0.007 


-0.082 


-0.143 


-0.039 


0.059 


-0.064 


-0.127 


-0.078 


0.033 


0.031 


0.067 


-0.438 


-0.532 


1.000 
















p^X. 


18 


0.024 


-0.052 


-0.045 


0.200 


0.098 


-0.118 


-0.119 


-0.198 


0.006 


-0.094 


-0.078 


-0.177 


-0.142 


0.028 


-0.445 


-0.235 


0.018 


1.000 














P^.. 


19 


0.022 


-0.002 


-0.033 


0.066 


0.557 


-0.077 


-0.104 


-0.617 


0.034 


0.003 


0.020 


-0.052 


0.259 


-0.377 


-0.237 


-0.020 


-0.089 


0.431 


1.000 












P% 


20 


-0.006 


-0.022 


-0.002 


0.084 


0.015 


-0.066 


-0.230 


-0.019 


0.011 


-0.013 


0.032 


-0.007 


-0.091 


0.001 


-0.017 


-0.001 


0.009 


0.060 


0.017 


1.000 










P%'" 


21 


-0.004 


0.113 


-0.071 


-0.088 


-0.002 


0.528 


0.069 


-0.000 


-0.063 


-0.201 


0.076 


0.023 


0.077 


0.005 


0.011 


0.080 


-0.085 


-0.048 


0.005 


-0.145 


1.000 








pir-^^ 


22 


0.002 


-0.109 


0.072 


0.080 


0.064 


-0.537 


-0.069 


-0.111 


0.061 


0.216 


-0.085 


-0.006 


-0.006 


-0.010 


-0.001 


-0.082 


0.076 


0.055 


0.061 


-0.514 


-0.671 


1.000 






p% 


23 


-0.004 


-0.005 


-0.025 


0.106 


0.040 


0.013 


-0.353 


-0.058 


-0.012 


-0.049 


0.043 


0.015 


-0.103 


-0.001 


-0.006 


0.015 


-0.011 


0.068 


0.038 


0.086 


0.027 


-0.028 


1.000 




P,o 


24 


0.012 


0.116 


-0.093 


-0.079 


0.020 


0.567 


0.156 


-0.036 


-0.082 


-0.232 


0.041 


-0.007 


0.075 


-0.002 


0.026 


0.071 


-0.097 


-0.045 


0.027 


-0.044 


0.304 


-0.294 


-0.064 


1.000 


OtherBB 

P." 


25 


-0.013 


-0.119 


0.103 


0.057 


0.022 


-0.598 


-0.067 


-0.043 


0.088 


0.257 


-0.058 


0.008 


-0.005 


-0.002 


-0.021 


-0.080 


0.096 


0.044 


0.021 


0.021 


-0.319 


0.319 


-0.472 


-0.781 1.000 



or 
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TABLE XVIII: Statistical (fit) correlation matrix of the B 
7v£u AB/Aq^ measurement for the four-mode fit. 



g^ range (GeV^) 0-4 4-8 8-12 


12-16 16-20 


>20 


0-4 1.000 0.191 0.050 


-0.005 0.068 


0.057 


4-8 1.000 0.089 


0.058 0.085 


0.011 


8-12 1.000 


0.197 0.127 


0.005 


12-16 


1.000 0.135 


-0.008 


16-20 


1.000 


0.032 


>20 




1.000 



TABLE XXII: Statistical (fit) correlation matrix of the B 
p£u AB/Aq^ measurement for the four-mode fit. 



g^ range (GeV^) 0-8 8-16 >16 



0-8 

8-16 

>16 



1.000 0.264 0.137 

1.000 0.189 

1.000 



TABLE XIX: Systematic correlation matrix of the B 
AB/Aq measurement for the four-mode fit. 



■kIv 



q^ range ( GeV^) 



0-4 



8-12 12-16 16-20 >20 



0-4 

4-8 

8-12 

12-16 

16-20 

>20 



1.000 0.521 0.705 0.394 

1.000 0.853 0.687 

1.000 0.652 

1.000 



■0.052 0.075 
0.605 0.478 
0.366 0.439 
0.637 0.367 
1.000 0.509 
1.000 



TABLE XX: Total correlation matrix of the B 
AB/Aq^ measurement for the four-mode fit. 



■kIv 



g^ range 


(GeV^) 


0-4 


4-8 


8-12 


12-16 


16-20 


>20 


0-4 




1.000 


0.337 0.401 


0.212 


0.015 


0.066 


4-8 






1.000 


0.430 


0.343 


0.272 


0.205 


8-12 








1.000 


0.443 


0.227 0.219 


12-16 










1.000 


0.350 


0.180 


16-20 












1.000 


0.221 


>20 














1.000 



TABLE XXI: Total correlation matrix of the B -i> -kIv 
AB/ A(f measurement for the four-mode fit after unfolding 
of the (f spectrum. 



g^ range ( GeV^) 0-4 4-8 8-12 12-16 


16-20 >20 


0-4 1.000 0.272 0.331 0.216 


-0.037 0.045 


4-8 1.000 0.390 0.273 


0.252 0.172 


8-12 1.000 0.475 


0.194 0.170 


12-16 1.000 


0.462 0.042 


16-20 


1.000 0.195 


>20 


1.000 



TABLE XXIII: Systematic correlation matrix of the B 
AB/Aq measurement for the four-mode fit. 



p£u 



g-" range (GeV^) 0-8 8-16 >16 



0-8 

8-16 

>16 



1.000 0.339 0.692 

1.000 0.296 

1.000 



TABLE XXIV: Total correlation matrix of the 
AB/Aq^ measurement for the four-mode fit. 

g^ range (GeV^) 0-8 8-16 >16 

0-8 1.000 0.307 0.532 

8-16 1.000 0.281 

>16 1.000 



B 



pe.v 



TABLE XXV: Total correlation matrix of the B -> plv 
AB/Aq measurement for the four-mode fit after unfolding 
of the q^ spectrum. 



g^ range (GeV^) 0-8 8-16 >16 



0-8 

8-16 

>16 



1.000 0.574 0.380 

1.000 0.389 

1.000 
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TABLE XXVI: Statistical (fit) covariance matrix of the B — )> 
7v£u AB/Aq^ measurement for the four-mode fit after unfold- 
ing of the q^ spectrum in units of 10~^^. 



g^ range (GeV^) 0-4 4-8 


8-12 


12-16 


16-20 


>20 


0-4 4.039 0.436 


-0.134 


-0.020 


-0.015 


0.116 


4-8 1.861 


0.135 


-0.027 


0.104 


0.018 


8-12 


1.462 


0.404 


0.110 


-0.018 


12-16 




1.720 


0.534 


-0.157 


16-20 






1.995 


0.014 


>20 








0.650 



TABLE XXVII: Systematic covariance matrix of the B — > nlv 
AB/Aq^ measurement for the four-mode fit after unfolding 
of the q^ spectrum in units of 10~^^. 



g^ range (GeV^) 0-4 4-8 


8-12 12-16 


16-20 >20 


0-4 4.425 0.888 


1.836 1.254 


-0.176 0.022 


4-8 0.931 


1.018 0.920 


0.653 0.289 


8-12 


1.666 1.244 


0.504 0.338 


12-16 


2.123 


1.091 0.245 


16-20 




1.228 0.359 


>20 




0.488 



TABLE XXVIII: Total covariance matrix of the B ^ ttIv 
AB/ A(f measurement for the four-mode fit after unfolding 
of the (f spectrum in units of 10~^''. 



g^ range ( GeV^) 0-4 4-8 8-12 12-16 


16-20 >20 


0-4 8.463 1.324 1.702 1.234 


-0.191 0.139 


4-8 2.792 1.152 0.894 


0.757 0.307 


8-12 3.129 1.648 


0.615 0.320 


12-16 3.843 


1.625 0.089 


16-20 


3.223 0.373 


>20 


1.138 



TABLE XXIX: Statistical (fit) covariance matrix of the B — >■ 
pHv AB/Aq^ measurement for the four-mode fit after unfold- 
ing of the q^ spectrum in units of 10~^^. 



g^ range (GeV^) 0-8 8-16 >16 



0-8 

8-16 

>16 



1.798 



0.536 
0.543 



■0.093 
0.066 
0.269 



TABLE XXX: Systematic covariance matrix of the B — > plv 



AB/Aq^ 
of the (f 



measurement for the four-mode fit after unfolding 
spectrum in units of 10^^^. 

q-" range (GeV^) 0-8 8-16 >16 



0-8 

8-16 

>16 



2.500 1.653 1.230 

2.837 0.968 

1.816 



TABLE XXXI: Total covariance matrix of the B ->■ plv 
AB/ Aq^ measurement for the four-mode fit after unfolding 



of the q spectrum in units of 10 



q-" range (GeV^) 0-8 8-16 >16 



0-8 

8-16 

>16 



4.298 



2.188 
3.381 



1.137 
1.034 
2.086 



